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Abstract 

A  software  package  was  written  which  will  analyze  an¬ 
nular  shaped  laser  wave-fronts.  The  polynomials  used  to 
estimate  the  wave-front  are  based  on  the  work  of  J.Y.  Wang 
and  D.E.  Silva  in  their  paper  "Wave-front  Interpretation 
with  Zernike  Polynomials."  This  involved  generating  a  set 
of  orthogonal  polynomials  from  the  Zernike  polynomials  by 
using  the  Gram-Schmidt  orthogonalization  process.  The  co¬ 
efficients  of  the  polynomials  are  determined  by  using  the 
orthogonality  of  the  polynomials,  instead  of  using  the  com¬ 
mon  least-squares  method.  The  coefficients  of  the  gener¬ 
ated  polynomials  arc  converted  to  Zernike  coefficients, 
and  both  sets  of  coefficients  are  presented  to  the  user. 

By  using  the  first  moments  of  the  wave-front's  posi¬ 
tion  in  the  collection  array,  the  software  is  able  to  de¬ 
fine  the  basic  parameters  of  the  wave-front.  These  param¬ 
eters  are:  center,  outside  radius,  and  the  obscuration  ra¬ 
tio  of  the  wave-front.  With  these  parameters,  the  soft¬ 
ware  computes  6,  then  11,  then  22  coefficients  to  show  the 
stability  of  the  coefficients.  With  well-defined  circular 
and  annular  wave-fronts,  the  program  was  consistantly  able 
to  compute  the  coefficients  with  an  RMS  error  of  less  than 
0.050  waves. 
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LASER  WAVE-FRONT  ANALYZER  SOFTWARE  IMPROVEMENT 


I.  Introduction 


The  Laser  Wave-front  Analyzer  (LWA)  at  the  Air  Force 
Weapons  Lab  (AFWL)  is  used  to  measure  the  phase  and  inten¬ 
sity  of  a  High  Energy  Laser  (HEL)  beam.  The  LWA  uses  the 
phase  and  intensity  data  to  calculate  the  coefficients  of 
Patterson's  polynomials  (Ref. 4).  The  coefficients  are  used 
to  determine  the  various  optical  aberrations  present  in  the 
laser;  the  aberrations  include:  piston,  tilt,  defocus,  and 
coma.  Once  the  aberrations  are  known,  they  can  be  used  to 
correct  the  laser's  optical  elements  to  maximize  the  far- 
field  intensity  of  the  laser. 

As  stated  earlier,  the  Patterson  polynomials  are  cur¬ 
rently  being  used  to  analyze  the  wave-front.  This  set  of 
polynomials  is  simple  and  is  not  the  best  set  of  polynomi¬ 
als  available.  A  better  set  would  be  the  Zernilte  polynom¬ 
ials  modified  for  an  annular  wave-front.  These  polynomials 
will  be  used  to  determine  the  coefficients  which  describe 
the  wave-front.  Since  the  polynomials  are  defined  over  an 
annular  wave-front,  the  coefficients  will  better  represent 
the  optical  aberrations. 


Problem 


This  thesis  involves  the  generation  of  a  FORTRAN  pro¬ 
gram  which  will  find  the  coefficients  to  both  Zcrnike  poly¬ 
nomials  and  a  set  of  polynomials  defined  over  an  annular  re¬ 
gion.  Since  most  IlEL's  have  annular  wave-fronts,  the  Zcrnike 
polynomials  must  be  modified  to  remain  orthogonal  over  the 
wave-front  (the  Zernike  polynomials  are  orthogonal  with 
circular  wave-fronts).  Orthogonality  of  the  polynomials 
will  refer  to  the  integral  over  the  entire  region  of  the 
product  of  two  polynomials.  When  the  polynomials  are  the 
same,  the  result  is  one  and  zero  when  they  are  not  the  same. 
The  coefficients  should  be  available  as  rapidly  as  possible 
after  the  collection  of  the  data  by  the  LWA. 

Assumptions 

Work  has  already  been  done  to  show  that  the  Zernike 
polynomials  are  not  orthogonal  when  part  of  the  beam  is  ob¬ 
scured  (annular)  (Ref.G);  therefore,  this  study  will  devel¬ 
op  the  software  required  to  evaluate  wave-fronts.  It  is  as¬ 
sumed  that  the  wave-fronts  are  annuli  whose  obscuration  ra¬ 
tio  ranges  from  0  to  almost  1.  The  obscuration  ratio  is  the 
ratio  of  the  inside  radius  to  the  outside  radius.  The  soft¬ 
ware  will  always  try  to  fit  the  wave-front  with  the  largest 
possible  annulus.  The  word  wave-front  will  mean  the  phase 
of  the  laser's  output,  at  any  instant  in  time.  Thus  with  an 
aberration-free  system  or  the  Zernike  coefficients  are  all 


p 


zero,  the  wave-front  will  be  planar. 

It  is  assumed  that  the  accuracy  of  the  results  will 
not  have  to  exceed  the  accuracy  of  a  16-bit  computer  (four 
significant  digits  to  the  right  of  the  decimal  point);  how¬ 
ever,  the  software  was  developed  on  a  60  bit  computer. 

Approach 

The  general  procedure  for  modifying  and  finding  the 
coefficients  to  the  Zernike  polynomials  has  been  presented 
in  Ref. 6  .  The  basic  procedure  consists  of  the  following 
steps:  (1)  remove  any  invalid  data  points,  (2)  find  the 
center  of  the  beam  and  determine  the  obscuration  ratio, 

(3)  using  the  Grarn-Schmidt  orthogonalization  method,  com¬ 
pute  the  modified  Zernike  polynomials,  (4)  solve  for  the 
aberration  coefficients,  and  (5)  print  the  results.  Each 
of  these  steps  will  be  presented  in  detail  later. 

Development 

The  basic  operation  of  a  LWA  and  the  polynomials  used 
to  analyze  wave-fronts  are  presented  in  Chapter  II.  The 
mechanical  and  software  operation  of  the  LWA  is  shown.  In 
Chapter  III,  the  conversion  of  Zernike  polynomials  to  a 
set  which  are  orthogonal  over  an  annulus  is  presented. 
Chapter  IV  deals  with  the  development  of  the  software. 

The  validation  of  the  software  is  presented  in  Chapter  V, 
and  the  conclusions  and  recommendations  are  presented  in 
Chapter  VI . 
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II .  Detailed  Analysis 


This  chapter  deals  with  the  basic  workings  of  a  Laser 
Wave-front  Analyzer.  Before  explaining  the  mathematics  re¬ 
quired  to  analyze  a  wave-front,  the  basic  operation  of  the 
LWA  will  be  discussed.  The  next  two  sections  present  an 
overview  of  some  of  the  methods  used  to  analyze  wave-fronts 
Included  in  the  sections  is  an  explanation  of  the  limita¬ 
tions  of  these  methods. 

LWA  Operation 

The  LWA  consists  of  two  basic  parts,  a  sliding  refer¬ 
ence  interferometer  to  collect  the  data  and  a  computer  to 
control  and  analyze  the  data.  The  inteferometer  is  design¬ 
ed  to  work  at  two  wavelength  bands:  9.1pm  to  10.7pm  and 
3.2pm  to  3.8pm  (Ref. 7:58).  The  LWA  is  capable  of  collect¬ 
ing  the  phase  gradient  and  intensity  of  the  beam  at  each 
point  of  a  32  by  32  array  100  times  a  second.  The  data  is 
collected  by  two  cryogenically  cooled  Hg:Cd:Te  detectors, 
one  for  the  X  axis  data  and  one  for  Y  axis  data.  The  phase 
map  is  assembled  later  under  software  control. 

The  actual  phase  measurement  is  done  by  collecting  the 
difference  in  phase  between  a  reference  beam  and  the  input 
beam.  The  input  beam  is  passed  through  the  optical  system 
shown  in  Figure  1,  and  is  focused  on  two  sets  of  apertures. 


-  V  S - 


*.* 


One  set  is  used  to  find  X  axis  data,  and  the  other  Y  axis 
data,  with  each  set  consisting  of  two  small  circular  aper¬ 
tures.  The  apertures  generate  Airy  diffraction  patterns 
and  their  proximity  to  each  other  generates  interference 
patterns  similar  to  Young's  double  slit  experiment.  The 
optical  system  focuses  the  interference  pattern  on  a  chop¬ 
ping  wheel,  which  is  used  to  modulate  the  two  beams.  The 
position  of  the  interference  pattern  is  proportional  to 
the  optical  path  difference  of  the  wave-front  at  the  two 
apertures.  Figure  2  shows  examples  of  parallel  and  tilted 
wave-fronts.  The  chopper  in  Figure  1  modulates  the  inter¬ 
ference  pattern  into  a  sine  wave.  Since  a  tilted  wave¬ 
front's  interference  pattern  is  shifted  off-axis,  its  mod¬ 
ulated  signal  will  have  a  corresponding  phase  shift.  In 
order  to  measure  the  phase  shift,  a  reference  beam  is  pass¬ 
ed  through  the  same  set  of  apertures  and  is  focused  on  its 
own  detector,  as  shown  in  Figure  1.  Since  the  modulated 
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reference  beam's  signal  is  constant  with  respect  to  time, 
the  phase  difference  between  it  and  the  input  beam's  sig¬ 
nal  is  used  to  represent  the  optical  path  difference.  The 
LWA  is  capable  of  detecting  a  phase  difference  as  small  as 
A/200  (Ref. 7:60).  The  phase  difference  is  called  delta 
phase.  The  delta  phase  data  is  collected  at  each  point  in 
the  array  by  using  two  orthogonal  pair  of  apertures,  re¬ 
sulting  in  X  and  Y  axis  delta  phase  maps.  At  the  same  time 
the  intensity  of  the  wave-front  is  measured  by  using  the  DC 
value  of  the  modulated  signal. 

Due  to  the  physical  layout  of  the  LWA,  the  X  and  Y 
delta  phase  is  not  collected  at  the  same  time.  The  X 
delta  phase  is  collected  starting  at  row  one  and  ending 
with  row  32.  The  Y  delta  phase  starts  at  column  24,  de¬ 
creasing  column-wise  to  one,  then  to  column  32  down  to  col¬ 
umn  25  (Ref. 2:6).  This  collection  method  is  alright  as 
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long  as  the  wave-front  does  not  change  shape  faster  than 
the  sample  period  of  0.01  seconds;  otherwise,  the  phase 
measured  at  one  pair  of  apertures  may  have  changed  by  the 
time  the  other  pair  of  apertures  is  at  the  same  point. 

The  delta  phase  and  intensity  data  is  digitized  such 
that  the  intensity  is  sent  as  an  eight-bit  word  and  the 
two  delta  phases  are  sent  as  eleven-bit  words,  with  one  bit 
used  to  signify  valid  data.  This  data  is  sent  as  two  16- 
bit  words  every  lOysecs  to  a  high  speed  storage  disk  to 
hold  the  data  until  it  can  be  processed.  The  next  step  in 
the  wave-front  analysis  is  the  reconstruction  of  the  wave- 
front.  This  is  discussed  in  the  next  section. 

LIYA  Software  Operation 

Once  the  data  has  been  collected  by  the  interferometer, 
the  computer  portion  of  the  LWA  must  first  reconstruct  the 
wave-front  from  the  delta  phase  data,  and  analyze  the  wave- 
front.  The  analysis  includes  isometric  and  isocontour  plots, 
far-field  intensity  distribution,  and  polynomial  coefficient 
fitting.  This  thesis  deals  only  with  the  coefficient  fitting. 

Before  presenting  the  current  polynomial  fitting  rou¬ 
tine,  the  method  used  to  reconstruct  the  wave-front  will  be 
presented.  The  basic  technique  is  to  use  a  least-squares 
technique  with  the  delta  phase.  Figure  3  shows  a  diagram  of 
a  basic  3  by  3  wave-front  collection  matrix.  Each 
point  on  the  array  corresponds  to  a  specific  value  of  the 
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Figure  3.  Simple  Wave-front  Collection  Array 
wave-front.  The  a^'s  and  b^'s  correspond  to  the  measur¬ 
ed  delta  phase  in  the  X  and  Y  direction  respectively. 

At  each  point  in  the  matrix  the  sum  of  the  square  of  the 
differences  is  computed.  For  example  at  the  point  W5 
the  difference  equation  is 

o5  =  (We-Ws-as)2  +  (W5-W.,-a2)2  +  (W8-W5-b,)2  +  (W5-W2-b2)2, 

Each  of  the  a  equations  are  differentiated  with  respect 
to  .  All  of  the  equations  are  presented  in  Appendix 

A.  For  this  example,  the  result  is  nine  equations  which 
can  be  solved  simultaneously.  As  a  reference,  one  of  the 
points  is  set  to  zero.  The  resulting  matrix  equation  is 

A  W  =  B  (2.2.1) 

where  W  is  a  column  matrix  on  the  actual  phase  values  of 
the  wave-front,  B  is  formed  from  the  measured  delta  phase, 
and  A  is  formed  from  the  difference  equations.  In  this 
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and  all  other  matrix  equations-;,  A  represents  a  n  by  n 
matrix,  A  represents  a  n  by  1  column  matrix,  and  A 
represents  a  scalar.  It  has  been  shown  by  others  (Ref. 5) 
that  A  in  Equation  (2.2.1)  is  Ilermitian,  irreducibly 
diagonally  dominant,  and  all  diagonal  entries  are  positive 
real  numbers. 

The  solution  is  slightly  more  involved  than  inverting 
A  and  multiplying  times  B  to  get  W  .  Some  of  the 
points  in  the  delta  phase  may  be  invalid  because  of  the  an¬ 
nular  region  of  the  wave-front.  The  software  computes  B 
by  trying  to  find  dominant  rows  and  columns  which  circum¬ 
scribe  the  annular  region.  Once  the  obscured  region  is 
"boxed"  in,  the  program  works  in  toward  the  inside  region 
and  outward  to  the  outer  radius  as  shown  in  Figure  4.  This 
is  done  to  reduce  the  possible  errors  introduced  by  edge 
effects . 

The  result  of  this  routine  is  a  reconstructed  wave- 
front  which  can  then  be  described  in  terms  of  various  poly¬ 
nomials.  The  software  used  now  utilizes  Patterson's  poly¬ 
nomials;  more  familiar  to  those  in  the  field  of  optics  are 
the  Zernike  polynomials.  The  next  two  sections  deal  with 
those  polynomials. 

Patterson  Polynomials 

The  Patterson  polynomials  are  a  set  of  polynomials 
that  were  developed  at  Perkin-Elmer  Corporation.  They  are 
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A  Central  Region  of 
Annulus 

B  Major  Axes  of  "Box" 

C  Direction  of  Scan 

D  Outside  Radius  of 
Annulus 

Figure  4.  Wave-front  Construction  from  Delta  Phase  Data 
basically  a  simpler  set  of  the  Zcrnike  polynomials.  Table 
I  shows  a  comparison  of  the  first  23  Patterson  and  Zernike 
polynomials.  The  most  important  point  about  the  Patterson 
polynomials  is  they  are  not  orthogonal  over  a  unit  circle; 
thus,  the  coefficients  of  the  polynomials  will  change  as 
more  or  fewer  of  the  different  polynomials  are  used  in  esti¬ 
mating  the  wave-front.  This  is  the  phenomenon  the  users  are 
faced  with  at  this  time.  As  they  increase  the  number  of  co¬ 
efficients  being  solved  for,  the  previously  found  coeffi¬ 
cients  change;  therefore,  they  do  not  know  when  the  coeffi¬ 
cients  are  the  actual  coefficients. 

Besides  using  polynomials  which  are  not  orthogonal,  the 
current  software  assumes  that  the  wave-front  will  fill  the 
32x32  collection  array,  and  that  it  is  also  centered  on  the 
same.  Thus  an  off-center,  small  beam  will  introduce  much 
more  error  than  the  non-orthogonality  of  the  Patterson  poly¬ 
nomials.  Figure  5  shows  this  case.  The  Patterson  polynomials 
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Table  I 


Polynomial  Comparison 

Patterson  (Ref  4)  Zernike  (Ref  8) 


P(l) 

1 

Z(l) 

1 

P(2  ) 

r  cos 

9 

Z(2) 

2  r  cos  0 

P(3) 

r  sin 

0 

Z(3) 

2  r  sin  6 

P(4) 

r2 

Z(4  ) 

S5  (2r2-l) 

P(5) 

r2cos 

20 

Z(5) 

/E  r2  sin  20 

P(G) 

r2  sin 

20 

Z(  6  ) 

/E  r2cos  20 

P(7) 

r3cos 

0 

Z(7) 

/5  (3r3-2r)sin0 

P(8) 

r3sin 

0 

Z(8) 

SE  (3r3-2r)cos0 

P(9) 

r3cos 

30 

Z(  9  ) 

/E  r3sin  30 

P(10) 

r3sin 

30 

Z(10) 

/g  r3cos  30 

P(ll) 

r1* 

Z(ll) 

/E  (6r“-6r2+l) 

P(  12  ) 

r4  cos 

20 

Z(  12  ) 

/10  (4r4 -3r2 )cos  20 

P(  13  ) 

r4  sin 

20 

Z(  13) 

/TO  (4r4 -3r2 )sin  20 

P(14) 

r1'  cos 

40 

Z(  14  ) 

/TO  r4cos  40 

P(15) 

r4  sin 

40 

Z(  15  ) 

/TO  r4sin  40 

P(1G) 

r5  cos 

0 

Z(16) 

/12  ( 10r5-12r3+3r)cos9 

P(17) 

r5  sin 

0 

Z(  17) 

/T2  ( 10r5 -12r3  +3r )sin0 

P(18) 

r5  cos 

30 

Z(  18) 

/12  (5r5~4r3)cos  30 

P(19) 

r5  sin 

30 

Z(19) 

/l2  (5rs-4r3 )sin  30 

P(  20 ) 

r5  cos 

50 

Z(20) 

/T 2  r5cos  59 

P(21) 

r5  s  in 

50 

Z  (  2 1 ) 

/T2  r5sin  50 

P(  22  ) 

r6 

Z(  22 ) 

/7  (20r6 -30r4  +12r2 -1) 
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Figure  5.  possible  Error  Condition  for  Input  Beam 


will  give  a  wave-front  which  appears  to  be  correct,  via  the 
RMS  error,  but  if  the  same  coefficients  were  used  to  recon¬ 
struct  the  actual  wave-front  with  the  proper  coordinates  for 
the  center  and  the  true  radius,  the  RMS  error  would  be  much 
higher. 

In  conclusion,  the  Patterson  polynomials  are  best  suit¬ 
ed  to  give  only  a  rough  approximation  of  the  wave-front  when 
it  fills  the  array  and  is  circular.  As  soon  as  the  wave- 
front  becomes  either  annular  or  does  not  properly  fill  the 
collection  array,  the  Patterson  polynomials  do  not  even 
closely  represent  the  aberrations  in  the  wave-front. 

Zernike  Polynomials 

The  Zernike  polynomials  are  the  classical  method  used 
in  the  field  of  optics  to  describe  the  aberrations  present 
in  a  wave-front  described  by  a  rotationally  symmetric  op¬ 
tical  system.  As  pointed  out  in  the  previous  section, 
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ito' 


these  polynomials  are  orthogonal  over  a  unit  circle.  Ap¬ 
pendix  B (p.105)  gives  the  reader  an  idea  of  what  the  first 
22  polynomials  look  like.  As  an  example,  Figure  6  shows 
the  results  when  the  first  five  coefficients  have  the  value 
of  0.1  waves,  where  one  wave  equals  one  wave-length  of  the 
wave- front . 

The  Zernike  polynomials  are  well  behaved  when  the  wave- 
front  is  circular.  But  when  the  wave-front  becomes  an  an¬ 
nulus,  the  Zernike  polynomials  are  no  longer  orthogonal 
over  the  annular  region.  Thus  they  are  not  well  suited  to 
the  analysis  of  most  HEL  beams.  Even  though  the  polynomials 
are  not  well  suited,  they  are  the  most  familiar  to  those  in 
optics;  therefore,  it  is  desirable  to  express  the  wave- 
front  in  terms  of  Zernike  coefficients.  The  next  chapter 
will  present  a  method  of  converting  the  coefficients  of  the 
newly  described  polynomials  to  the  coefficients  of  Zernike 
polynomials . 

Summary 

In  this  chapter  the  basic  operation  of  the  LWA  has  been 
presented,  along  with  the  way  the  wave-front  is  analyzed. 

The  previous  discussion  has  shown  that  neither  the  Zernike 
or  the  Patterson  polynomials  are  well  suited  to  define  the 
output  of  most  HEL's.  The  possible  sources  of  error  that  re¬ 
sult  from  the  sampling  and  analysis  have  also  been  presented. 
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The  next  chapter  will  present  a  method  of  generating  a 
set  of  polynomials  which  are  orthogonal  over  an  annulus. 
Not  only  will  they  bo  orthogonal  over  any-sized  annulus, 
but  the  coefficients  to  the  new  set  of  polynomials  can  be 
converted  to  Zernike  coefficients. 


15 


III. 


Thcor  o  tical  Developmen  t 


As  pointed  out  in  the  last  chapter,  the  Zernike  poly¬ 
nomials  are  not  orthogonal  over  a  region  defined  by  an  an¬ 
nulus.  This  chapter  deals  with  the  theory  behind  and  the 
development  of  a  set  of  polynomials  which  are  orthogonal 
over  an  annular  region  with  the  basic  Gram-Schmidt  method 
of  orthogonalization  being  presented  first.  This  will  then 
be  applied  to  the  Zernike  polynomials,  which  are  defined  in 
more  detail,  to  generate  the  new  set  of  polynomials  orthog¬ 
onal  over  an  annular  domain  .  Finally,  the  new  polynomials 
are  tested  to  verify  their  orthogonality.  This  chapter  is 
based  on  the  work  of  Y/ang  and  Silva  (Ref  8). 

Gram-Schmidt  0  r  t  ho  go  n  a  1  i  z  a.  t  io  n 

The  Gram-Schmidt  orthogonalization  process  is  a  method 
in  which  a  set  of  linearly  independent  vectors  is  combined 
to  form  a  new  set  of  vectors  which  arc  orthogonal.  In  par¬ 
ticular,  this  process  will  be  applied  to  the  Zernike  poly¬ 
nomials.  This  will  bo  done  in  the  next  section;  this  sec¬ 
tion  deals  with  the  specifics  of  the  Gram-Schmidt  orthogon¬ 
alization  process. 

The  basic  theorem  for  the  Gram-Schmidt  orthogonaliza¬ 
tion  process  is  to: 

let  the  vectors  u i , . . . , u  ,  form  a  set  of  linearly  in¬ 
dependent  set  in  a  vector'  space  V  with  inner  product 


(•,•)  .  Then  one  can  construct  an  orthogonal  set  of 

vectors  Xi,...,x  such  that  for  each  i  with  1-i-s 
the  set  of  vectors  xlt...,x.  spans  precisely  the  same 
subspace  as  does  the  set  of  vectors  (Kef. 3: 

138).  1 

The  basic  process  involved  to  generate  an  orthogonal  set  is 

to  take  the  first  vector  ux  and  to  normalize  it,  or 

vj  =  ui  ,  xi  =  - - - r  (3.1.1) 

( v  i  ,  v  i  ) ;’ 

where  (v  V  )  is  the  inner  product  of  v.  and  v.  .  To 
i  ’  .)  i  J 

span  the  same  vector  space,  the  second  vector  u2  has  a 
multiple  of  x  subtracted  from  it,  or 

v2  =  u2  -  HiXi  (3.1.2) 

where  oti  is  chosen  such  that  X]  and  V2  are  orthogonal. 
To  find  a  value  for  oti  ,  the  vector  xi  is  multiplied 
times  both  sides  of  Equation  (3.1.2)  and  the  inner  product 
is  then  taken  to  satisfy  the  orthogonality  requirement  of 
xi  and  V2  giving 

(X],vz)  =  ( x i , u 2 )  -  a i (x i , x ! )  (3.1.3) 

When  two  vectors  are  orthogonal,  their  inner  product  is 
zero.  Another  important  point  is  that  the  inner  product  of 
a  normalized  vector  with  itself  is  one.  Since  Equation  (3. 
1.1)  is  the  process  used  to  generate  a  normalized  vector, 
and  Xj  and  vz  must  be  orthogonal  to  satisfy  the  theorem, 
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Equation  (3.1.3)  becomes 


a i  =  ( x i , u  2 )  (3.1.4) 

Putting  this  result  into  Equation  (3.1.2)  gives 

V2  =  U2  -  (x ,  ,  u2  )x  j  ,x2  = - — — (3.1.5) 

(v2 , v2 )a 

The  third  vector  v3  follows  the  same  process,  thus 

v3  =  u3  -  a2x2  -  BiXi  (3.1.6) 

where  a2  and  Bi  are  found  such  that  v3  is  orthogonal 
to  Xj  and  x2 

As  an  exampi e,  just  & )  will  be  solved  for,  but  the 
process  is  the  same  for  a2  .  Thus 

(x1(v3)  =  (xi,u3)  -  a2(xj,x2)  -  3i(xlfxi)  (3.1.7) 

where  again  (x1(u3)  must  be  zero  to  satisfy  the  theorem 
and  (xi,x2)  is  zero  because  they  are  two  previously  de¬ 
fined  orthogonal  vectors.  Equation  (3.1.7)  becomes 

3i  -  ( x i , u 3 )  (3.1.8) 

and  in  the  same  manner 

ct 2  =  ( x 2  , u 3 )  (3.1.9) 

Using  Equations  (3. 1.8, 9)  in  Equation  (3.1.6)  gives 


This  process  is  continued  until  all  of  the  desired  x  vec¬ 
tors  are  produced  (Ref. 3: 138). 


*  • 


**  * 


( 


As  point  out  earlier,  the  above  process  is  the  basic 
means  to  produce  a  set  of  orthogonal  vectors.  It  has  been 
found  that  this  process  can  be  numerically  inaccurate  (Ref. 
3:141),  and  thus  a  modified  Gram-Schmidt  process  has  been 
developed.  In  the  original  method,  x2  was  computed  from 
ui  and  u->,...,u  were  left  alone,  with  each  new  x. 
vector  using  only  and  not  affecting  the  remaining  vec¬ 
tors  u.  u,  .  In  the  modified  Gram-Schmidt  process, 

besides  computing  x2  from  u2  ,  (x2,u.)x2  is  subt act¬ 

ed  from  ui  to  produce  a  new  ui  with  i  ranging  from 
two  to  S  to  make  u^  orthogonal  to  x2  (Ref. 3: 142). 
Since  u|  is  the  same  as  Equation  (3.1.5),  u^  is  normal¬ 
ized  to  produce  x2  .  As  before  this  vector  is  used  such 
that  (x2,uj)x2  is  subtracted  from  u^  to  produce  u? 
with  i  ranging  from  three  to  S  to  make  u?  orthogonal 
to  x2  as  well  as  x2  (Ref. 3: 142).  This  is  continued  un¬ 
til  xg  has  been  generated. 

The  following  example  shows  how  the  set  {x1,x2fx3} 
is  generated  from  {u2,u2,u,}  .  As  before 


(Ui ,Ui )“ 

but  the  remaining  vectors  are  modified  such  that 


(3.1.11) 


i 

U2  -  U2 


(Xj ,u2 )x i 


(3.1.12) 
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* 


(Xi ,u3 )xj 


(3.1.13) 


U  3  =  U  3  “ 


Again  from  Equation  (3.1.5) 


i 

u2 

~7~i  nr 

(u2  ,u2) 


which  is  used  to  modify  u3  giving 


2 


i 

=  u3 


( x  2 ,u3 )x2 


and 


(3.1.14) 


(3.1.15) 


2 

x3  =  - (3.1.16) 

(us-u3)2 

Notice  the  difference  between  Equations  (3.1.16)  and  (3.1. 
10).  In  Equation  (3.1.10)  the  inner  product  is  taken  with 
an  unmodified  x2  and  u3  .  Thus  the  modified  Gram- 
Schmidt  process  gives  more  accurate  results  over  that  of 
the  basic  process.  It  is  this  modified  Gram-Schmidt  pro¬ 
cess  which  is  used  to  generate  the  new  polynomials.  Before 
presenting  the  modification  of  the  Zernike  polynomials,  the 
basic  formulation  of  the  Zernike  polynomials  is  presented. 

Zernike  Polynomials 

The  basic  characteristics  of  the  Zernike  polynomials 
were  discussed  earlier,  whereas  this  section  deals  with 
their  basic  construction.  Table  I  showed  22  of  the  Zernike 
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-C.. 


Zernlke  Radial  Polynomials 


polynomials.  It  can  be  seen  that  each  polynomial  consists 
of  the  product  of  a  radial  and  an  angular  term.  Since  the 
new  set  of  polynomials  deals  with  annular  regions,  the  poly¬ 
nomials  will  be  independent  of  any  angular  functions;  there¬ 
fore,  it  is  desirable  to  separate  these  two  terms. 

The  means  of  deriving  the  radial  terms  of  the  Zernike 
polynomials  has  already  been  done  by  others  (Ref . 1 : sec . 9. 2 ) 
The  basic  process  is  described  in  Appendix  C.  Table  II 
shows  the  radial  polynomials  and  all  of  the  corresponding 
Zernike  polynomials.  Each  radial  polynomial  was  produced 
by 


R-(r)  = 


(n-m)/2 

£ 

s=0 


( -l)s(n-s) !  rn~2s 


s!  [(n+m)/2  -  sj!^(n-in)/2  -  sjl 


(3.2.1) 


where  m  and  n  are  radial  degree  and  azmuthal  frequency 
respectively.  The  following  conditions  must  also  be  met 
with  m  and  n  ,  such  that  m-n  and  n  minus  m  is  even 
The  Zernike  polynomials  are  generated  by  the  following 
set  of  equations: 


Z 

even  j 

Zodd  j 


£2(n+l)ja  R™(r)  cos  m0 
[2(n+l)]  ^  R™(r)  sin  m0 


m^O " 

► 


Z.  =  [(n+1 

3  1 


)p  R1’ 

'  i  n 


(r) 


m=0 


(3.2.2) 


22 


where  j  is  the  mode-ordering  number  (Ref .8 : 1510) .  In 
other  words,  if  the  second  Zernike  polynomial  is  desired, 
j  is  set  to  two.  Besides  the  conditions  on  m  and  n 
mentioned  earlier,  the  following  properties  of  the  Zernike 
polynomials  exist: 

(1)  The  polynomials  are  invariant  with  respect  to  ro¬ 
tation  about  the  center  of  the  unit  circle. 


(2)  R™(r)  =  R"m(r) 


(3)  /R™(r)R^,(r)rdr  = 


(3.2.3) 


(3.2.4) 


(4) 


J  d2rW(r)Zj(r)Zj.(r)  = 


(3.2.5) 


where 


W(r)  = 


-c 


for [ r | -1 
for | r I <0 


and 


/!  i=J 
Sij  =■) 

U  ifj 


The  first  ten  Zernike  polynomials  have  other  names 
which  are  familiar  to  those  in  the  field  of  optics.  Table 
III  lists  these  names  with  the  corresponding  Zernike  poly¬ 
nomial  .  With  the  definition  of  the  Zernike  polynomials 
which  are  orthogonal  over  an  annular  region;  this  is  the 
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Table  III 

Aberration  and 

Corresponding 

Zernike  Polynomial 

Z(l) 

Piston 

2(2  ) 

X  Tilt 

Z(3) 

Y  Tilt 

2(4 ) 

Defocus 

Z(  5  ) 

0°  Astigmatism 

Z(6  ) 

45°  Astigmatism 

2(7) 

Y  Coma 

2(8) 

X  Coma 

Z(9) 

Y  Clover 

Z(10) 

X  Clover 

topic  of  the  next  section. 

Generation  of  Orthogonal  Polynomials 

In  the  previous  section,  the  equations  which  define 
the  Zernike  polynomials  were  presented,  along  with  the 
specific  properties  of  the  Zernike  polynomials.  In  par¬ 
ticular,  Equation  (3.2.5)  deals  with  the  orthogonality  of 
the  Zernike  polynomials  and  the  region  of  space  in  which 
they  arc  valid.  In  the  case  of  an  annulus,  the  weighting 
function  W(r)  can  no  longer  have  the  domain  of  -1  to 
1  but  must  have  a  domain  of  6-|rj-l  ,  where  B  is  the 

obscuration  ratio. 
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Since  an  annulus  is  independent  of  any  angular  func¬ 
tion,  the  orthogonal ization  procedure  need  only  deal  with 


the  radial  portion  of  the  Zorn  ike  polynomials.  Thus  the 
modified  Crain-Schmi  dt  procedure  will  only  be  performed  on 
the  radial  portion  of  the  Zorn  ike  polynomial.  Again  re¬ 
ferring  to  the  aspect  of  angular  independence,  each  new 
polynomial  N™(r)  will  only  depend  on  those  radial  poly¬ 
nomials  with  the  same  azmuthial  frequency  m  ,  or  Nm(r) 

m 

will  depend  on  NTm(r)  ,  Nm,2(r)  »•••*  Nn_2(r)  * 

o 

Starting  with  R0(r)  the  modified  Gram-Schmidt  pro¬ 
cedure  yields 


o  K0(r) 

N0(r)  =  — - - - - -r 

(H0(r),R0(r))“ 

where 

1 

(R°(r) ,R°(r) )  -  f  Rj(r)R°(r)rdr 

*4 


(3.3.1) 


(3.3.2) 


which  comes  from  the  definition  of  the  inner  product  of  two 
continuous  functions.  Since  R°(r)  has  no  lower  order  ra¬ 
dial  terms  (n^O)  ,  this  is  the  final  step  for  this  new 

polynomial.  This  is  also  true  for  every  new  polynomial 
where  tji  equals  n  ,  since  one  of  the  conditions  on  m 
and  n  with  the  Zernike  radial  polynomials  required  that 
n  must  be  greater  or  equal  to  m  .  Thus  in  general 
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Nn(r) 
n v  ' 


Rn  ( r ) 


(  Kn  (  r  )  Rn  (  r  )  ) :! 

v  n  '  >  n 


(3.3.3) 


where  again  in  general 


( Rn  ( r  >  ’  Rk ( 1  )  >  =  f  Hn(r)Rk(r)rdr 


(3.3.4) 


and  where  rn  ,  n  ,  j  and  k  are  integers  and  obey  the 

rules  set  forth  in  the  previous  section.  Since  I'k(r)  is 

already  covered  in  Equation  (3.3.3),  the  next  polynomial  is 

N°(r)  .  In  this  case  the  vector  space  is  R°(r)  and 

2  o 

li°(r)  .  N°(r)  is  found  the  same  as  that  in  Equation 

2  o 

(3.3.1)  but  R°(r)  must  be  modified  to 
2 

R°'(r)  =  R2°(r)-(N;(r),R20(r))N0°(r)  (3.3.5) 


and  therefore 


N2(r)  = 


K  '(r) 


(R2°'(r),R2°'(r))a 


(3.3.6) 


2 

Again  Equation  (3.3.3)  covers  the  case  of  N2(r)  .  Since 

1  o 

N3(r)  is  very  similar  to  N2(r)  ,  the  next  new  polynomial 

to  be  developed  will  be  N°(r)  .  This  polynomial  will  use 

0  0  0 

as  its  vector  space  H0(r)  ,  R2(r)  ,  and  R4(r)  .  The 

process  starts  just  the  same  as  with  N2(r)  with  Equations 

(3. 3. 1,5,6).  The  next  step  is  to  modify  R°(r)  with  N°(r) 
o 

and  then  with  N2(r)  ,  thus 
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(3.3.7) 


R°  ''(I')  =  Rj(r)-(Nj(r),Rj(r))Nj(r) 


when  N0(r)  was  computed  in  Equation  (3.3.1)  and 


R°~(r)  =  u“'(r)-(N2°(r),l{iC/(r))N"(r) 


(3.3.8) 


when  N°(r)  was  computed  in  Equation  (3.3.6). 
2 


Finally 


N.  (r)  = 


R„"(r) 


(11  ~(r),R  ~(r)) 


(3.3.9) 


This  process  is  continued  until  all  of  the  desired 


N™(r)'s  have  been  found.  Obviously  xf  22  polynomials  are 
to  be  used  in  the  wave-front  analysis,  the  process  of  com¬ 
puting  all  of  the  integrals  would  be  time-consuming.  It 
is  possible  to  reduce  this  time  factor  by  realizing  that 
once  6  is  found,  all  of  the  inner  product  integrals  can 
be  done  at  once.  The  next  section  deals  with  the  genera¬ 
tion  of  these  integration  terms. 


Inner  Product  Coefficients 

As  mentioned  in  the  previous  section,  the  fastest  way 
to  compute  the  new  radial  polynomials  would  be  to  compute 
all  of  the  inner  product  integrals  at  once.  This  process 
would  have  to  be  done  in  an  orderly  fashion  since  each  in¬ 
tegral  is  dependent  on  the  previous  integral.  For  example 
in  Equation  (3.3.9),  the  inner  product  terms  are  dependent 
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1  wv* 


-*  -  on  the  results  of  the  inner  product  .integrals  in  Equations 

(3.3. 1,6).  It  would  appear  that  each  integral  would  have 
to  be  computed  before  doing  the  next  one.  Mathematically 
this  is  true,  but  since  this  system  of  dependent  equations 
will  be  put  into  a  program,  the  problem  can  be  solved  in 
a  different  way.  First  a  set  of  variables  which  will  ease 
the  computation  of  solving  for  the  new  radial  polynomials 
will  be  derived. 

In  Equation  (3.3.1),  if  the  inner  product  in  the  de¬ 
nominator  is  changed  to 

v° 

Yoo  =  (R°(r),R°(r))a  (3.4.1) 

o  o 

and  by  substituting  Equation  (3.4.1)  in  Equation  (3.3.1) 
one  gets 

o  R°(r) 

N  (r)  »  _2 -  (3.4.2) 

Y#„ 

0  o 

m  m 

where  ynn  is  the  inner  product  of  ETn(r)  and  Rj(r)  for 

j=m,m+2 , . . . , n-2  ,  and  the  square  root  of  the  inner  product 

of  R™(r)  anc*  ^j(r)  wllon  n=j  .  To  clarify  this,  the 

radial  polynomial  N°(r)  will  be  recomputed  making  use  of 

4 

the  y  terms. 

o  o  o 

As  before,  N  (r)  uses  the  vector  space  {ll0(r) ,  Rj  (r) , 

o 

Rl|(r)}  .  Thus  using  Equation  (3.4.2)  as  a  starting  point, 

the  remaining  two  vectors  must  be  modified  yielding 


R 0  "  ( r )  -  R0(r)-(N0(r),R0(r))N°(r) 
2  2  0  2  0 


(3.4.3) 


and 


R 0  "  (  r  )  •■=  R°(r)-(N0(r),R0(r))N°(r) 

4  4  0  4  0 


(3.4.4) 


The  inner  product  terms  of  Equations  (3.4. 3,4)  become 


Y°  =  (N° (r) ,R° (r) ) 
2  0  0  2 


(3.4.5) 


and 


Y°  =  (N° (r ) , R° (r )  ) . 

4  0  0  4 


(3.4.6) 


Using  Equations  (3.4. 5,6)  in  Equations  (3. 4. 3, 4)  yields 


RB/(r)  =  R°(r)-Y  N°(r) 

2  2  2  0  U 


(3.4.7) 


and 


R°  (r )  =  R°(r)-Y°  N°(r) 

4  4  4  0  0 


(3.4.8) 


Continuing  with  the  modified  Gram-Schmidt  procedure  gives 


N°(r)  = 


R°'(r) 

2 


(R°^(r),R°'(r))“ 
2  2 


(3.4.9) 


with  the  inner  product  term  becoming 


Y°  =  (H0-(r),R°'(r))“ 
2  2  2  2 


(3.4.10) 
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With  N°(r)  computed,  N°(r)  can  lie  found  by 


li  (r)  -  Ii('"(r)-(Nc(r),R0"(r))N0(r) 

It  4  2  4  2 


(3.4.11) 


with 


4  2 


(N°(r),U°'(r)) 

2  4 


(3.4.12) 


and 


N°(r) 

4 


R°-(r) 

4 


(ll°-(r),R0-(r)): 

4  4 


(3.4.13) 


The  final  inner  product  term  is 


y°  =  (R0(r),R0~(r))* 

4  4  4  4 


(3.4.14) 


The  y  terms  can  be  generalized  with 


i 

■l 


m  _  l  Mm  .  . nm .  .  , 

Y  .  =  [  N  ,.(r)R  (r)rdr 
nj  )  m+ j  v  7  nv 


(3.4.15) 


where  j=0 , 2 , 4 , . . . , n-2  and 


m 

^iin 


r  1 

/< 

n 


(N™*(r))2rdr 


(3.4.16) 


where 


».m ^ .  .  _  ..m.  .  m  ,.m  ,  . 

N  (r)  =  R  (r)-Y  „N  0(r)- 
n  v  '  n  'nn -2  n-2v  ' 


-Ym  Nm(r)  (3.4.17) 
'  nm  m  v  ' 
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From  Equation:;  (3.4.15,10),  the  entire  set  of  y's  can  be 
computed . 

It  appear;.:  that  all  the  y  terms  do  is  simplify  the 
equations,  but  if  Equation  (3.4.13)  is  expanded  in  terms 
°i  ( r ) ' s  from  Equations  (3.4.1-12),  the  resulting  equa¬ 
tion  is 


N 0  ( r  ) 


1 

~B — 
4  4 


>\(r)-  ( i’)-  ~a—(.  E 0  ( r )  -  -J-°-R°(r)) 


0  0 


2  2 


o  a 


(3.4.18) 


Solving  Equation  (3.4.18)  for  R°(r)  in  terms  of  y's  and 

4 

N™(r)'s  gives 

R°(r)  ~y 0  N°(r)+y°  N°(r)-iy°  H°(r)  (3.4.19) 

4  444  422  400 


Therefore  each  Zernike  radial  polynomial  can  be  expressed 
in  terms  of  the  new  radial  polynomials  and  vice  versa.  Tn 
general  Equation  (3.4.19)  becomes 


nm,  \  _  m  ..no  .  ,  m 

R  ( r )  =  y  N  (r)+y 
n^  '  nm  nr  nm+2 


( r ) + . . .+ym  Nm(r) 
nn  n  ' 


(3.2.20) 


In  a  previous  section  the  basic,  requirements  of  the 
Zernike  polynomials  wore  expressed.  Those  requirements, 
being  slightly  modified,  should  still  be  met.  Thus 


Nn(r)  “  Nnn,(r> 


(3.4.21) 


> 


Mm,  ...in  .  .  .  (Snn 

N  (r)N  ^(r)rdr  =  rn — rrr 
n  n  7  2(n+.1) 


(3.4.22) 
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/  d  r  r  W  (  r )  jN  Z  .  ( r  )  N  Z  j  ^  (  r )  =  6..  , 


(3.4.23) 


whore 

fl/v  B-  |  r  |  >-l. 

w(2)  =< 

(O  otherwise 

and  NZ.(r)  is  the  newly  defined  Zorn ike  polynomial,  using 
the  N™(r)'s  as  the  radial  components.  Equation  (3.4.21) 
is  valid  since  N™(r)  is  made  up  of  R™(r)'s  which  meet 
this  requirement  as  shown  in  Equation  (3.2.3).  The  next 
section  deals  with  the  meeting  of  the  remaining  require¬ 
ments. 


Orthogonality  Tests 

As  pointed  out  in  the  previous  section,  it  is  desir¬ 
able  for  the  new  radial  polynomials  N™(r)  to  have  the 
same  qualities  as  the  Zernike  polynomials  R™(r)  •  In 

particular  the  orthogonality  of  the  radial  polynomials  are 
shown  in  Equations  (3.2.3)  and  (3.4.22)  for  the  Zernike 
and  the  new  radial  polynomials  respectively. 

As  a  test  of  the  new  radial  terms,  the  first  test  will 
be  with  the  polynomials  where  m=n  ,  thus  Equation  (3.4.22) 


becomes 


1 

f N”(r )n” ^(r )rdr 

*B 


L  Ann " 

2  (n+.l ) 


(3.5.1) 


32 


Using  Equation  (3.4.17),  the  two  integrals  cancel,  giving 
the  proper  results;  thus,  Equation  (3.5.7)  is  valid  for 
all  values  of  m-n 

The  final  test  is  to  see  if  Equation  (3.5.8)  still 


holds  true  when  it  is  changed  to 


(3.5.11) 


where  n^.j  .  Since  j  can  not  equal  n  the  Kronecker 

delta  (Snj)  will  be  zero;  therefore,  Equation  (3.5.11) 

should  be  zero.  This  integral  is  very  difficult  to  verify 

using  the  general  case,  since  each  new  radial  polynomial 

Nm(r)  and  Nin(r)  would  have  to  be  expanded  similar  to 
n  j 

Equation  (3.4.  ID).  Since  in  ,  n  ,  and  j  are  arbi¬ 
trary,  the  series  becomes  too  unwieldy  to  determine  the 
point  at  which  errors  may  have  been  introduced;  therefore, 
several  test  cases  were  developed  using  specific  values 
for  m  ,  n  ,  and  j  .  Since  the  process  is  still  very 
tedious,  only  a  simple  case  and  the  final  results  will  be 
presented. 

As  a  simple  test  case 

1 

/N°(r)N°(r)rdr  2  o  (3.5.12) 

2  0 

e 


will  be  tested.  Each  of  the  polynomials  can  be  expanded  to 


make  Equation  (3.5.12)  become 


j  -£-(K“(r)-  --fiRo(r)) 

g  2  2  0  0 


-4-li°(r)  rdr  2  0 
Y  o 
o  o 


(3.5.13) 


From  Table  II  and  simplifying,  this  becomes 


/*  Y  9 

/  2r2 -( 1+  -f-0-)  rdr  =  0 

el  Y°°- 


(3.5.14) 


Using  Equations  (3.4.15,16)  and  integrating  yields  £(32-B4) 
which  does  not  equal  zero.  Since  the  equation  used  to  gen¬ 
erate  y°  has  already  been  verified,  the  only  error  can 
o  o 

come  from  the  generation  of  y°  .  As  a  first  guess  Equa- 

2  0 

tion  (3.4.15)  is  changed  to 


Yjj  j  “  2(n+l) 


m.  .(r)Rm(r)rdr 
m+jv  n v  ’ 


(3.5.15) 


but  this  gives  32-BM  when  used  in  Equation  (3.5.12).  The 
next  change  made  worked,  making 


Y"j  =  2 ( m+ 1 ) 


f  *m*  ■ 

J  ni+J 


(r)R“‘(r)rdr 


(3.5.16) 


This  equation  was  found  invalid  when  the  test  case  using 
N°(r)  and  N°(r)  was  used.  The  final  form  of  the  equa- 

>t  2 

tion  is 
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This  equation  was  verified  using  several  test  cases,  which 
are  too  tedious  to  present.  With  Equations  (3.4.20),  (3.5. 
7),  and  (3.5.17),  it  is  possible  to  generate  all  of  the  new 
radial  polynomials;  thereby,  using  N™(r)  P"*-ace  °f  R™(r) 

in  the  Zernike  polynomials,  a  new  set  of  orthogonal  polynom¬ 
ials  can  be  generated.  Figures  7  to  11  show  the  inner  pro¬ 
duct  coefficients  with  respect  to  the  obscuration  ratio. 

With  the  set  of  orthogonal  polynomials,  an  annular  wave- 
front  can  now  be  more  accurately  expressed.  The  method  on 
analyzing  the  wave- front  will  be  presented  next. 

Wa vc- f ron t  Analysis 

As  presented  by  others  (Kef . l:Soc.9. 1,2) ,  wave-fronts 
can  be  represented  as  the  summation  of  the  Zernike  polyno — 
mials  times  their  respective  coefficients,  or 

oo 

ajZj(r)  =  <f>  (3.6.1) 

j  =  l 

where  <j>  is  the  measured  wave-front,  and  a.  is  the  coeffi- 

J 

cient  to  the  Zernike  polynomial  (Z.(r)).  Usually  in  optics 

J 

one  is  concerned  with  the  lower  order  Zernike  polynomials; 
therefore,  Equation  (3.6.1)  becomes, 


o 
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Figure  9.  Inner  Product  Coefficients 
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Inner  Product  Coefficients 


Inner  Product  Coefficients 


(3.6.2) 


N 


£ 


x.7.  .(r) 

.)  J 


*  4' 


j  =  l 


As  pointed  out  in  the  last  chapter,  t)>  is  not  a  continuous 
function,  but  is  sampled  at  discrete  points.  Thus  Equation 
(3.6.2)  can  bo  further  modified  to 


N 


£ 


a  .  Z  .  (  r  .  ) 
J  J  3 


4>(ri) 


(3.6.3) 


where  i=l,2,...,M  ,  and  M  is  the  number  of  sampled 

points.  This  equation  can  bo  written  in  the  matrix  form 


:(r.)a  «  (KiV) 


(3.6.4) 


where  as  before  a  is  a  N  by  1  column  matrix,  Z(r.) 

ti  1 


is  a  1  by  N  row  matrix,  and  (j>(r^)  is  a  scalar. 


The  most  common  method  of  solving  for  the  coefficients 
is  to  use  a  least-squares  approach.  This  involves  minimiz¬ 
ing  the  sum  of  the  squares  of  the  difference  between  the 
estimated  and  the  measured  wave-fronts,  or 


M 


A  = 


£ 


1=1 


£  V^'d-V^d 

j  =  l 


(3.6.5) 


To  find  each  coefficient,  the  derivative  of  A  is  taken 


with  respect  to  each  a^  and  sot  to  zero.  Thus 
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*1 


0 


(3.6.6) 


wr 


M 

N 

“ 

9  A  \p 

9a  r  zL/ 

k  i=l 

2i5h(?i) 

Z  .  (r  .  )a  r  .  ) 

L-j  j  i  j  i 

_ 

or 


M  N 

Z  (r.  )  Z.(r.  )a  . 

Z_^  k  i  Zw  J  1  J 

j  =  l 


i-1 


i-1 


( 3 .  G  .  7 ) 


This  equation  can  be  written  in  matrix  form  as 


T  T 

Z  Z  a  =  Z  (J> 


(3.6.S) 


and  following  matrix  inversion  and  multiplication  one  gets 

(3.G.9) 


T  - 1  T 
a  =  (ZXZ)  Z  4 


11  11  IT 


The  problem  with  this  method  for  finding  the  coefficients 
T 

is  the  matrix  Z  Z  can  be  numerically  unstable;  therefore, 

ii  it 

the  correct  solution  of  the  equation  is  difficult,  but  not 
impossible  to  find  (lief .  8  : 1514 ) . 

The  method  used  in  this  thesis  to  find  the  coefficients 
of  the  polynomials  uses  the  orthogonality  of  the  polynomials, 
As  before;,  the  least-squares  method  is  used  giving 


M 

i=l 


N 


7  b.NZ.(r.)-<J>(r.) 

A-j  J  J  i  i 

U:1 


(3.6.10) 


where  b.  is  the  new  polynomial's  coefficient.  Equation 
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. 


•  * 


(3.4.23)  can  be  approximated  in  discrete  form  as 


I 


NZ  (r.)NZ  (r.)  =  6 
.W  U  I  Vv  l'  uv 


1-1 


(3.6.11) 


Taking  the  derivative  of  Equation  (3.G.10)  and  setting  it 
equal  to  zero  gives 


3  A 

3b. 

k 


M 


E 

i-1 


N 


j“i 


=  0  (3.6.12) 


Since  NZ  (r  )  is  a  constant  with  respect  to  the  summation 
over  j  ,  it  can  be  put  inside  the  second  summation.  Us¬ 
ing  Equation  (3.6.11),  the  sum  from  j~l  to  N  is  non-zero 
only  when  j  equals  k  ;  therefore,  Equation  (3.6.12)  be¬ 
comes 


b 


k 


1 

M 


M 


ENZk(fi)+(?i) 

i-1 


(3.6.13) 


With  this  equation,  the  coefficients  can  be  found  directly 
with  no  matrix  inversion.  Notice  the  same  is  true  for  the 
Zorn  ike  polynomials  when  the  region  is  circular. 

From  this  section  and  the  previous  one  it  has  been 
shown  that  it  is  possible  to  construct  a  set  of  new  poly¬ 
nomials  which  are  orthogonal  over  an  annular  region.  These 
new  polynomials  can  be  used  to  find  the  coefficients  to 
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describe  an  annular  wnve-l'ront  without  employing  the  stand¬ 
ard  matrix  inversion  method.  'J’his  makes  the  programing  much 
more  reliable  and  faster.  Since  the  now  polynomials  are  de¬ 
rived  from  the  Zorn  ike  polynomials,  the  coefficients  can  be 
transposed  into  the  corresponding  Zcrnikc  coefficients. 

This  process  will  be  explained  in  the  next  chapter  which 
deals  with  the  process  of  developing  software  to  utilize  the 
results  of  this  chapter. 
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IV.  Program  Dovel opmont 


This  chapter  explains  the  development  of  software 
which  performs  the  analysis  of  annular  cave-front  using 
the  theory  presented  in  the  last  chapter.  The  software 
consists  of  two  major  sections,  the  definition  of  the  re¬ 
gion  of  analysis  and  the  actual  analysis  of  the  wave-front. 
This  chapter  presents  the  main  program  and  each  of  the  maj¬ 
or  subroutines,  with  the  subroutines  being  presented  in  the 
order  in  which  they  are  called.  Each  section  will  present 
tiie  basic  algorithm,  a  basic  flowchart,  and  a  description 
of  the  algorithm  as  necessary. 


Main  Progr am 


The  main  program  is  the  controlling  program  of  most  of 
the  subroutines.  It  does  the  actual  analysis  of  the  wave- 
front  and  prints  the  results  for  the  user.  The  basic  al¬ 
gorithm  for  this  program  is: 

1.  Read  in  first  frame  of  data  and  call  PAPER.  FAPER 
returns  the  computed  radius,  center,  and  obscuration 
ratio.  The  operator  is  asked  if  he  wants  to  change 
any  parameter  and  the  number  of  frames  to  bo  analyzed. 

2.  If  not  the  first  time  through  the  routine,  then 
read  in  a  new  frame  of  data. 

3.  Define  region  of  analysis  as  given  from  parameters 
computed  or  entered  in  step  1  by  calling  CONTUR. 

4.  Do  steps  f>  to  11  three  times:  once  to  compute  6  co¬ 
efficients,  then  to  compute  11  coefficients,  then  to 
compute  22  coefficients. 
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5.  Zero  out  all  data  arrays  used  in  computation. 

G.  Using  the  obscuration  ratio,  compute  the  inner 
product  coefficients  by  calling  GAMSUJ1. 

7.  Compute  X  and  Y  starting  positions  and  increments. 

8.  Using  the  region  defined  in  step  7  and  the  param¬ 
eters  from  step  -1,  compute  the  new  polynomial  coeffi¬ 
cients  . 

9.  Convert  coefficients  from  step  8  to  Zernikc  coeffi¬ 
cients. 

10.  Generate  two  estimated  wave- fronts  from  the  Zernike 
and  new  polynomial  coefficients. 

11.  Compute  RMS  error  and  print  results. 

12.  If  not  done  with  all  frames  do  to  step  2. 

13.  Stop. 

This  algorithm  is  presented  in  flowchart  form  in  Figure  12. 
Since  each  subroutine  is  presented  later,  only  those  steps 
which  do  not  rely  on  major  subroutines  are  presented  here. 

The  first  step  in  the  algorithm  uses  FAPER  to  define 
the  wave-front.  After  this  is  done,  the  user  is  asked  if 
he  wants  to  see  the  defined  region  and  the  computed  param¬ 
eter.  These  parameters  include  the  center  of  the  annular 
region,  the  outside  radius,  and  the  obscuration  ratio.  The 
user  is  then  asked  if  he  wants  to  change  any  parameters  if 
he  knows  the  actual  values.  This  request  is  made  to  correct 
any  possible  errors  that  FAPER  may  have  made.  The  reasons 
for  error  in  FAPER  will  be  presented  later.  The  program 
will  also  ask  the  user  how  many  frames  he  wants  to  have 


Figure  12.  Basic  Flowchart  of  Main  Program 
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analyzed.  With  Llii.s  data,  the  program  is  ready  to  start 
the  actual  analysis. 

Alter  the  .first,  frame  has  been  analyzed,  a  new  frame 
is  read.  This  frame  of  data,  is  defined  into  regions  of  no 
data,  valid  data,  and  the  center  of  the  annulus  by  the 
parameters  determined  from  the  first  frame.  This  is  done 
by  t!ie  subroutine  CONTUK.  Once  the;  regions  of  the  wave- 
front  are  defined,  the  program  starts  to  analyze  the  wave- 
front.  Each  wave-front  has  6,  then  11,  then  22,  coeffi¬ 
cients  computed  to  verify  that  the  number  of  coefficients 
does  not  affect  the  values  of  the  coefficients.  Before 
each  new  set  of  coefficients  is  computed,  the  data  arrays 
are  zeroed. 

Using  the  obscuration  ratio  parameter,  the  subroutine 
GAMSU13  is  called.  This  subroutine  computes  the  inner  pro¬ 
duct  coefficients  needed  to  convert  Zorn ike  radial  poly¬ 
nomials  to  the  new  set  of  radial  polynomials.  The  result¬ 
ing  lower  triangular  matrix  GAMMA  is  inverted  and  placed 
into  INGAMMA  for  future  use.  Each  array  is  used  in  the  con¬ 
version  process  depending  on  which  radial  polynomial  is  be¬ 
ing  converted  to  the  other.  With  these  arrays,  the  program 
can  begin  to  analyze  the  phase  data. 

Since  the  phase  array  is  not  always  filled,  the  range 
of  X  and  Y  depend  on  the  size  of  the  radius.  In  other 
words,  if  the  radius  is  only  ten  units  in  a  32  by  32  unit 
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Figure  13.  Possible  Cartesian  Coordinates  of  Phase  Array 
array,  the  position  of  the  wave-front  will  change  the  X 
and  Y  values  of  the  four  corners.  This  can  be  seen  in 
Figure  13. 

Using  the  radius  and  center,  the  program  must  deter¬ 
mine  the  X  and  Y  increment  and  the  initial  values  XO 
and  YO  .  The  analysis  starts  in  quadrant  II  of  the  coor¬ 
dinate  system,  and  scans  from  -X  to  X  ,  Y  to  -Y  end¬ 
ing  in  quadrant  IV  .  XO  and  YO  are  found,  such  that 


XO  = 


Radius 


(4.1.1) 


Y-l _ 

Radius 


(4.1.2) 


where  X  and  Y  are  the  center  coordinates  of  the  wave- 
front  given  in  terms  of  increments  in  the  32  by  32  array. 


Thu-;  i  ho  pn  i  n  t  (1,1)  in  the  array  corresponds  to  the  point 
(  -XO ,  YO )  in  (:t  i- 1  <  s  i  an  oaordinal  i's.  The  .i  ne  rernon  t  val  lies 
for  X  and  i  are  simply  foiuid  by  Ink  in;;  t  lie  reciprocal 
of  th1'  radius  n  !  the  \va  \  e-  iron  t. . 

The  lire.  I  step  in  the  analyai  a  is  to  compute  the  co¬ 
efficients  to  the  new  set  of  polynomials.  As  pointed  out 
in  the  last  chapter,  each  new  radial  polynomial  is  made  up 
of  tv  linear  combination  of  Zernike  radial  polynomials.  This 
can  be  expressed  in  matrix  form  as 

INZ  -  R  (4.1.3) 

where  y_  is  a  square ,  lower  tr iangular  matrix,  NZ  is  the 
new  set  of  radial  polynomials,  and  R  is  the  Zernike  radial 
pol  ynomittl  s ,  It  was  pointed  out  that  either  of  the  radial 
polynomials  are  used  with  the  angular  dependence  to  construct, 
the  tot  a. I  polynomial;  therefore,  the  program  computes  the  an¬ 
gular  data  and  multiplier,  it  times  either  radial  polynomial 
to  generate  Live  total  polynomial.  Thus  to  compute  the  new 
polynomial,  the  following  matrix  equation  is  used, 

NZ  =  Cy~ *  R  (4.1.4) 

where  C  is  a  diagonal  array  of  angular  constants.  As  an 
example  the  equation  can  be  expressed  as 
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where  the  first  four  new  polynomials  are  desired. 

To  actually  compute  the  coefficients,  the  program  will 
scan  the  phase  array.  As  soon  as  the  radial  value  on  the 
point  in  the  scan  is  less  than  one,  /x2 +y2 -1  ,  the  Zernike 

radial  polynomials  R  and  the  angular  constants  C  are  com¬ 
puted  and  the  matrix  multiplication  indicated  in  Equation 
(4.1.4)  is  performed.  Using  Equation  (3.6.13),  the  coeffi¬ 
cient  for  each  polynomial  is  computed.  Once  the  entire  phase 
array  is  scanned,  the  coefficients  are  found  by  dividing  the 
summation  by  the  total  number  of  points  used  in  determing 
the  coefficient.  This  process  can  be  expressed  as 

M 

bk'sE  (4.1.5) 

i=l 

where,  as  before,  b^  is  the  coefficient  to  the  new  poly¬ 
nomial  NZ^(r^)  ,  and  M  is  the  number  of  valid  data  points. 
A  valid  data  point  is  one  which  is  inside  the  annulus  or  unit 
circle,  and  is  a  phase  measurement.  Thus,  using  Equation 
(4.1.4)  and  (4.1.5),  the  new  polynomial  coefficients  are  de¬ 
termined.  The  next  step  is  to  convert  these  coefficients  to 
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Zernike  polynomial  coefficients. 

The  Zorn  ike  coefficients  can  be  easily  found  by  first 
converting  the  new  polynomial  coefficients  to  coefficients 
of  the  new  radial  polynomials.  Next  the  radial  coefficients 
are  converted  to  Zorn  ike  radial  coefficients  and  finally  t.o 
Zernike  polynomial  coefficients.  This  process  can  be  ex¬ 
pressed  as 

a  CyC-)b  (4.1.6) 

where  b  is  the  new  polynomial  coefficients,  and  a^  is 
the  Zernike  polynomial  coefficients.  As  before,  the  process 
is  done  be  scanning  through  the  array  and  calculating  the 
values  of  C  and  C  1  at  each  valid  point.  Each  coeffi¬ 
cient  is  kept  as  a  running  number  of  valid  points  at  the  end 
of  the  scan.  Equation  (4,1.6)  can  bo  expressed  as 

M 

‘I'lE  (4.1.7) 

i=l 

With  the  two  sets  of  coefficients,  the  next  step  is  to  com¬ 
pute  the  estimated  wave-fronts. 

Reconstructing  the  estimated  wave-front  is  basically 
the  opposite  of  finding  the  coefficients.  The  only  differ¬ 
ence  is  that  the  two  wave-fronts  can  be  found  at  the  same 
time.  The  matrix  equations  are 

<b  =  a  C  R  (4.1.8) 

/j  — ~ 
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and 


‘f'nZ  “  llC>  (4.1.9) 

where  d>  and  6  are  the  wave- fronts  found  with  Zernike 
z  nz 

and  new  polynomial  coefficients  respect ively .  The  wave- 
fronts  are  formed  by  scanning  through  the  32  by  32  array  and 
at  each  valid  point,  the  matrix  Equations  (4. 1.8,9)  are  eval¬ 
uated.  Tlius  after  the  scan,  there  are  two  estimated  wave- 
fronts  . 

Once  the  estimated  wave-fronts  are  constructed,  the  Root 
Mean  Square  (RMS)  error  is  calculated.  This  error  is  found 
by  calling  the  subroutine  RMSERR.  Its  operation  is  explain¬ 
ed  later.  Once  the  RMS  error  is  found,  the  main  program  will 
print  out  the  coefficients  of  both  polynomials  and  the  RMS 
error  from  their  respective  estimated  wave-fronts.  The  pro¬ 
gram  will  then  either  return  to  compute  the  next  set  of  co¬ 
efficients,  or  start  a  new  frame  of  data. 

In  summary,  the  main  program  does  the  actual  analysis 
of  the  wave-front.  It  also  asks  for  any  input  parameters 
and  calls  various  subroutines  to  prep a  re  the  frame  for  an¬ 
alysis.  Appendix  U  contains  detailed  flowcharts  of  this 
routine  and  all  the  subroutines.  The  next  section  deals 
with  the  subroutine  which  computes  the  inner  product  coeffi¬ 


cients. 


(.ambit. 

Tit"  subruul  ine  GAMSl'B  is  a  foul  jmt  which  cumpu  t  cs  the 

IT) 

inner  product  co<*  1  r  i  c  i  out  s ,  y  .  This  routine*  has  throe* 

nj 

pa  ram*  *1  ors  in  its  cal  line,  statement:  the  array  in  hold  the 
inner  product  coefficients,  the  obscuration  ratio,  and  the 
number  of  coefficients. 

As  pointed  out  in  the  last  chapter,  the  actual  inner 
product  terms  are  dependent  on  each  other.  To  make  the 
program  run  faster  and  more  accurately,  each  of  the  integrals 
in  liquations  (3,5.7,17)  were  done  by  hand.  This  results  in 
a  series  of  dependent  equations  with  a  dependency  on  lower 
order  terms  and  the  obscuration  ratio,  p  .  Since  the  in¬ 
tegrations  where  done  by  hand,  each  integral  does  not  have 
to  bo  done  by  some  numerical  method,  i.e.  Simpson's  rule*  or 
the  trapezoidal  method.  This  greatly  reduces  the  time  since 
only  one  pass  through  the  routine  is  required  instead  of 
integrating  each  term.  The  error  is  reduced  by  not  approx¬ 
imating  a  term  which  is  used  in  the  approximation  of  another 
term. 

Since  this  subroutine  is  just  a  series  of  equations, 
the  algorithm  will  not  be  presented.  GAMSUD's  flowchart  is 
presented  in  appendix  D,  and.  the  entire  matrix  it  generates 
is  in  appendix  E.  Before  GAMSUB  can  run,  the  obscuration 
ratio  must  be-  found.  It  is  found  in  the  subroutine  PAPER 
which  is  presented  in  the  next  section. 
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The  .subroutine  TAPER  return.'.;  the  center,  outside  ra- 


F.Ai-i'li 


d  i  us 

,  and 

the  obsciii'; 

it  ion  ratio.  It  also 

defines  the  value 

for 

no  da 

la  present. 

Tills  subroutine  is 

the  most  important 

otto 

s  i  nci' 

an}  errors 

can  affect  where  the 

main  program 

"thinks"  the  wave -front  is  positioned.  The  main  purpose  of 
this  subroutine  is  to  define  the  type  of  wave-front,  i.e. 
annular  or  circular,  and  then  find  the  parameters. 

Since  most  of  the  analysis  will  be  with  annular  wave- 
fronts,  I-APTR  will  assume  an  annular  wave-front  is  present. 
The  algorithm  for  this  subroutine  is; 

1.  Select  the  value  to  be  used  as  no  data  (BAD)  as  the 
value  which  is  present  at  any  three  corners  of  the  col¬ 
lection  array. 

2.  If  the  point  (10,10)  in  the  array  contains  a  value 
of  BAD,  go  to  step  5. 

3.  If  not,  starting  at  (16,16),  spiral  outward  until 
a  point  with  a  value  of  BAD  is  round  or  the  radial 
distance  from  (16,16)  is  greater  than  5.75. 

4.  When  a  BAD  point  is  found,  set  now  center  to  coor¬ 
dinates  of  this  point.  If  the  distance  was  greater 
than  5.75,  then  keep  the  center  at  (16,16). 

5.  By  calling  EDGJi,  find  the  vertical  and  horizontal 
radii  of  the  inside'  radius  of  the  annulus.  If  EDGE 
returns  an  error  condition,  the  wave- front  is  non-annu- 
lar;  therefore  go  to  step  6.  Repeat  until  the  radii 
stabilize  or  ten  iterations  have  been  done. 

6.  From  the  center  found  in  step  4,  find  the  smallest 
outside  radius. 

7.  If  the  region  was  defined  as  non-nnnular,  define 
the  wave-front  by  calling  CO.NTUR  ,  and  find  the  largest 
outside  radius  which  will  tit.  inside  the  wave-front. 
Pass  this  value  and  the  center  coordinates  to  VALID  to 
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I'iud  tin.;  Floating-point  vuIucjj  of  11k?  radius  and  cun- 
t  or. 

8.  TF  the  region  if;  annular,  u.so  C0NTUR  to  define  the 
region.  VAhll>  Ls  called  to  lint!  the  smallest  inside 
radius  tor  the  annul  in;  which  will,  circumscribe  the1  HM) 
data.  The  inside  radius-  and  the  smallest.  outside  radi¬ 
us  is  used  to  do Tine  the  largest  possible  outside  radius. 

9.  Fill  the  character  array  with  final  symbols. 

10.  Return  to  calling;  routine. 

The  basic  flowchart  for  this  subroutine  is  shown  in  Figure 
14  . 

The  first  two  steps  of  the  algorithm  are  very  straight 
forward.  The  determining  of  the  value  of  no  data  is  crucial 
to  the  entire  program's  operation.  This  and  other  subroutines 
use  this  value  as  a  test  to  determine;  whether  to  use  a  point 
in  the  analysis  of  the  wave- front.  It  is  also  used  to  de¬ 
fine;  the  central  region  of  the  annulus.  The  next  step  is  to 
find,  and  if  possible,  define  the  inner  radius  of  the  annu¬ 
lus. 

When  the  LWA  is  run,  the  input  wave-front  is  usually  cen¬ 
tered  on  the  collection  array;  therefore,  this  subroutine  will 
start  looking  at  l ’no  center  For  the  annulus.  If  a  valid  data 
point  is  present  at  the  center  (115,1.0),  then  the  subroutine 
will  spiral  outward  to  try  to  find  the  central  region .  Fig¬ 
ure  15  shows  the  shape  and  direction  of  the  spiral  at  the 
largest  extent  of  the  search. 

If  the  data  at  any  point  along  the  search  is  invalid, 
assume-  it  has  Found  the  edge  of  the  inside 
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the  program  will 


Figure  15.  Spiral  Search  Pattern  of  PAPER 


region  of  the  annulus.  Invalid  means  that  the  point  has  a 
value  of  BAD.  The  maximum  radius  of  5,75  for  the  search 
was  chosen  as  an  optimum  value  from  several  test  cases. 

The  entire  search  process  is  based  on  the  assumption  that 
the  operator  of  the  EV.’A  has  aligned  the  system  to  fill  the 
array.  If  an  invalid  point  is  found  in  the  search,  this 
point  is  tentatively  set.  as  the  new  center  of  the  annulus. 

The  next  stop  is  to  find  the  approximate  center  of 
the  annular  region.  Assuming-  the  wave- 1  rent  is  annular, 
the  subrout  inc  EDGE  is  called  to  find  the  horizontal  and 
vertical  radii  from  l he  center  in  all  four  directions. 
Using  the  averages  1,1  t  he  !.«'<>  hor  i  son  I  a  1  and  two  vertical 
radii,  a  new  renter  is  computed  and  EDGE  is  called  again. 
This  process  is  repented  until  ei  liter  the  radii  do  not 

GO 


change  or  Lon  iterations  arc'  performed.  Tho  limit  of  ton 
iterations  is  to  prevent  the  possibility  of  the  integer 
center  oscillating  around  the  actual  center.  If  EDGE  runs 
into  the  limits  of  the  array,  then  EAPER  will  assume  the 
spiral  seal'd)  found  the  outside  edge  of  a  circular  wave- 
front.  If  this  happened,  the  center  is  reset  to  (16, 16). 

At  this  point  the  program  "knows"  either  the  approximate 
location  of  the  center  of  the  annulus  or  that  the  wave- 
front  is  non-annular. 

In  either  type  of  wave-front,  the  next  step  is  to  find 
the  smallest  horizontal  and  vertical,  radii  from  the  center. 
With  this  radius,  in  the  case  of  an  annulus,  VALID  is 
called  lo  find  the  center  and  tho  smallest  radius  which 
covers  tlu'  inside  region  of  the  annulus.  When  the  wave- 
front  is  non-annular,  the  radius  is  used  to  find  the  larg¬ 
est  radius.  VALID  is  again  called  to  look  for  the  center 
and  the  largest  circle  which  will  fit  in  tho  region. 

Jf  the  wave'- front  .is  non-annular,  then  PAPER  returns 
the  results  o  •'  VALID  lo  the  calling  routine.  Otherwise 
PAPER  wiJJ  search  for  tin'  outside  radius  of  the  annulus. 
This  is  done  by  using  the  floating-point  values  of  the 
center  found  by  VALID.  This  radius  and  the  inside  radius 
are  used  to  compute'  the  obscuration  ratio,  BETA.  BETA,  the 
center,  and  the  outside  radius  are  then  returned  to  the 
calling  routine. 
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I u  summary,  f’APER  is  the*  controlling  routine  to  define 
the  typo  and  parameters  of  the  wave-front.  This  routine  is 
called  to  show  the  operator  what  it  "thinks"  the  wave-front 
looks  like.  With  this  Information,  the  operator  ca.n  accept 
or  reject  its  results;  thereby  reducing;  the  likelihood  of 
error.  The  next  section  deals  with  the  subroutine  VALID. 


VALID 


The  subroutine  VALID' s  main  purpose  is  to  compute  the 
floating-point  value  of  the  center  and  the  radius  of  the 
wave-front.  As  pointed  out  in  the  previous  section,  VALID 
has  two  possible  modes  of  operation:  find  the  center  and 
outside  radius  of  a  non-annular  wave-front,  or  find  the  cen 
ter  and  inside  radius  of  an  annular  wave-front.  The  basic 


algorithm  is: 

1.  Using  the  flag  passed  to  it,  go  to  step  3  if  a 
non-annular  wave-front . 

2.  Define  the  region  of  search  to  the  center  ±  the 
estimated  radius  plus  one.  Go  to  step  4. 


3.  Define  the  search  area  as  the  entire  array. 


4.  Search  the  defined  region  for  either  valid  or  in¬ 
valid  points  depending  on  whether  the  wave-front  is 
non-annular  or  annular. 


5. 

and 

if 


Keep  a.  running  total  of  the 
the  values  of  the  vertical 
looking  for  the  center. 


number  of  points  found 
and  horizontal  indices 


G.  Divide  the  sums  of  the  vertical  and  horizontal  in¬ 
dices  by  the  total  number  of  points,  to  get  the  true 
center  of  the  region. 


7.  Using  the  new  center  and  the  estimated  radius, 
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repent  the  search.  If  annular  region,  count  the  num¬ 
ber  of  invalid  points,  decreasing  the  radius  of  search 
by  0.2b.  As  soon  as  the  new  count  changes  from  the 
previous  count,  use  the  last  value  of  the  radius  as 
the  smallest  radius.  If  non-annular,  increase  the 
radius  until  an  invalid  point  is  found  and  use  previous 
radius  as  the  largest  radius. 

8.  In  non-annular,  return  the  new  center  and  the  radius 
to  the  calling  routine. 

9.  If  annular,  fill  character  array  SEU  with  the  sym¬ 
bol  for  the  central  region  of  the  annulus  and  return 
the  calling  routine. 


A  flowchart  o  ('  this  routine  is  shown  in  Figure  16. 

As  indicated  in  the  algorithm,  and  the  previous  sec¬ 
tion,  VALID  is  passed  the  integer  estimate  of  the  center 
and  radius  from  the  calling  routine.  It  uses  these  values 
to  try  to  find  the  floating-point  values  of  the  same.  To 
avoid  confusion,  the  algorithm  will  be  split  into  the  two 
cases  of  annular  and  non-annular  analysis.  In  the  case  of 
an  annular  wave-front,  VALID  will  use  the  integer  radius 
to  circumscribe  the  inside  region  of  the  annulus  with  the 
search  region.  On  the  first  pass  through,  VALID  will  com¬ 
pute  the  first  moments  of  the  X  and  Y  values  of  the  co¬ 
ordinates  of  the  invalid  data  points.  The  moments  are  found 
by  using  the  following  equations 


XC  =» 


i-1 


v  .X 


(4.4.1) 
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and 


E 


v  .  Y 

1 


where 


v  . 


x 


{1  when  the  data  is  invalid 
0  otherwise 


(4.4.2) 


NT  is  the  number  of  invalid  points  in  the  search  area,  and 
N  is  the  number  of  points  in  the  entire  search  area.  XC 
and  YC  in  these  equations  are  the  floating-point  values 
of  the  center.  Once  the  center  is  found,  the  search  is  re¬ 
peated  and  the  smallest  radius  is  found  which  will  include 
all  of  the  invalid  data  points  (NT).  The  radius  is  changed 
in  steps  of  0.25.  This  decrement  was  selected  because  the 
use  of  a  floating-point  radius  and  center  in  a  discrete 
array  gives  a  reasonable  decrement  without  sacrificing 
accuracy.  Once  the  radius  is  found,  VALID  fills  the  char¬ 
acter  array  SLR  with  the  symbol  to  represent  the  cen¬ 

ter  of  the  annulus.  VALID  then  returns  its  results  to  the 
calling  routine. 

In  the  case  of  a  non-annular  r  gion,  VALID  defines  the 
search  area  as  the  entire  data  array.  It  then  follows  the 
same  procedure  as  before.  VALID  first  finds  the  center  by 
using  the  first  moments  as  in  Equations  (4.4. 1,2)  but 

v  Jl  when  the  data  is  valid 

i  )0  otherwise 
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Figure  17.  Error  Caused  by  VALID* s  Means  of  Finding  Center 

Using  the  computed  center  and  the  estimated  outside  radius, 
VALID  will  search  the  array.  If  an  invalid  point  is  found 
in  the  radius  of  search,  the  radius  :i  s  decreased  by  0.25. 
This  repeated  until,  no  invalid  points  are  found;  thus,  the 
smallest  radius  is  found  which  circumscribes  the  valid  data. 
The  radius  and  center  are  then  returned  to  the  calling 
routine . 

Since  VALID  uses  the  first  moment  to  find  the  center, 
errors  can  be  introduced.  Figure  17  shows  one  such  example 
where  the  input  beam  is  a  semicircle.  This  is  a  case  where 
user  intervention  is  important.  The  user  can  reset  the  cen¬ 
ter  and  radius  parameters  to  best  fit  the  wave-front.  VALID 
works  best  with  wave-fronts  which  are  circularly 
symmetric.  The  further  the'  wave-front  is  from  being  symme¬ 
tric.,  the  greater  the  possibility  of  error. 


G6 


The  subroutine  VALID  is  used  to  find  the  floating¬ 


point  values  of  the  center  and  radius,  from  integer  guesses. 
It  works  best  with  wave-fronts  which  are  symmetric.  VALID 
can  give  erroneous  data;  therefore,  the  main  program  will 
check  with  the  user  before  continuing. 

EDGE 

The  subroutine  EDGE  is  called  by  PAPER  to  estimate  the 
center  of  the  annulus.  EDGE  has  two  modes  of  operation: 
vertical  or  horizontal  scan.  The  mode  is  selected  by  a 
flag  which  is  passed  to  EDGE.  The  algorithm  for  EDGE  is: 

1.  Starting  at  coordinates  givc-n  by  calling  routine, 
search  left  or  down  until  a  valid  point  is  found. 

Save  the  distance. 

2.  Starting  at  the  same  coordinates ,  search  right  or 
up  until  a  valid  point  is  found.  Save  (die  distance. 

3.  If  the  search  in  any  direction  hits  the  outside 
boundary  of  the  data  array,  then  set  error  flag  and 
return.  Otherwise  return  left/right  or  up/down  dis¬ 
tances  to  calling  routine. 

The  basic  flowchart  for  tins  algorithm  appears  in  Figure  IS. 

As  pointed  out  in  the  section  on  PAPER,  EDGE  is  called 
at  the  most  ton  times  to  try  to  find  the  best  guess  for  the 
center  of  the  annulus.  The  distance  returned  is  also  used 
in  the  determination  of  the  inside  radius.  The  next  sections 
deal  with  subroutines  which  are  used  in  the  analysis  of  the 


wave-front . 


ZKAD 


The  subroutine  7, HAD  is  used  to  generate  the  values  of 
the  Zerniko  radial  polynomials  at  a  specific  X  ,  Y  , 
coordinate.  liquations  (4. 1.4,5)  show  that  to  find  the  co¬ 
efficients  to  new  polynomials,  the  new  polynomials  must  be 
computed  at  each  X  ,  Y  point.  ZKAD  is  a  very  simple 
subroutine,  as  shown  in  the  algorithm.  The  algorithm  is: 

1.  Using  the  X , Y  values  and  the  number  of  terms 
needed,  compute  the  Zerniko  radial  polynomial  terms 
Rm(X,Y). 

2.  Fill  a  column  matrix  with  the  desired  number  of 
terms . 

ZKAD  uses  the  radial  polynomials  presented  in  Table  II. 

As  an  example  when  six  coefficients  are  needed,  ZKAD  will 

produce  a  column  matrix  consisting  of  {K° (X , Y) , K1 (X ,Y ) , 

0  1 

K1(X,Y),R0(X,Y),K:'(a,Y),K--(X,Y)}  .  The  repetition  of 

12  2  2 

some  of  the  terms  conies  from  Equat  ion  (3.2.2).  The  index 
of  the  matrix  K  corresponds  to  j  in  Equation  (3.2.2). 
The  entire  matrix  generated  by  ZKAD  is  presented  in  appen¬ 
dix  E.  Thus  ZKAD  is  a  simple  but  necessary  routine'  for  the 
generation  of  either  Zerniko  or  new  polynomials. 

zanc; 

The  subroutine  ZANG  is  similar  to  ZKAD.  ZANG  computes 
the  angular  dependence  of  the  Zerniko  polynomials  for  a 
specific  X,Y  coordinate.  Instead  o  T  producing  a  column 
matrix,  ZANG  makes  a  square  diagonal  matrix.  The  matrix 

<>9 


is  .shown  in  Append  ix  K .  The  index  into  the  army,  ( j  ,  j  ) 
where  j  \1  to  the  number  of  coeffi  ci  ent:s ,  i.s  I  he  mode- 
order  ,j  in  Equation  (3.2.2).  ZANG  simply  computes  the 
constant  and  s i n. /cosine  portion  of  the  Zernike  polynom¬ 
ials.  ZANG  and  ZKAD  are  separate  from  each  other  to  allow 
easy  conversions  of  new  polynomial  coefficients  to  Zernike 
coefficients  as  shown  in  Equation  (4.1.6). 

CONTUK 

The  subroutine  CONTUK  is  used  to  fill  a  character  array 
with  symbols.  CONTUK  uses  the  center  coordinates  and  the 
inside  radius  of  the  annulus  to  fill  the  array  SER  with: 

when  the  data  point,  is  valid,  when  no  data  is 

present.,  and  for  the  center  of  the  annulus. 

CONTUK  is  a.  simple  program  which  scans  the  phase  data, 
putting  in  the  proper  symbol.  When  the  scan  gets  inside  the 
central  region  of  the  annulus,  the  is  used  regardless 

of  the  value  of  the  data  point. 

CONTUK  is  the'  last  subroutine  which  is  used  in  the  an¬ 
alysis  of  the  data.  The  next  sections  deal  with  either  data 
manipulation  or  the  output  of  results. 


KVSEKK 

The  .subroutine  KMSEKR  finds  the  normalized  RMS  error 
between  two  arrays.  The  basic  equation  used  is 
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ERROR 


▼  H 


where  A.  is  the  input  wave- I' runt  array,  13.  is  the 
1  i 

estimated  wave-front ,  and  N  i s  the  total  number  of  points 
used  in  computing  the  estimated  wave-front.  To  find  the 
number  of  points,  RMSERR  uses  the  "A"  in  the  character 
array.  This  symboi  was  placed  in  the  array  by  F APISH  to 
telJ  the  remaining  routines  whore  valid  data  for  analysis 
can  be  found.  Thus  whenever  a  "A"  is  found,  RMSERR  then 
adds  another  term  to  the  summation.  The  normalization  al¬ 
lows  a  better  comparison  between  obscured  wave-fronts  and 
n on -  obs cu v o d  v,  a  v o -  f rents. 


Minor _ Subs out i nos 

The  minor  subroutines  include  SERPRNT,  INVERT,  MULT, 
ARPR1,  and  ARPR2.  The  subroucino  SERPRNT  is  called  when¬ 
ever  the  character  array  SER  is  printed.  SERPRNT  also 
prints  out  the  parameters  used  in  the  analysis  of  the  wave- 
front.  The  parameters  are:  the  center,  radius,  and  obscur¬ 
ation  ratio. 

INVERT  and  MULT  are  two  matrix  routines.  INVERT  in¬ 
verts  a  lower  triangular,  square  matrix.  It  is  INVERT 
which  inverts  the  inner  product  coefficient  array  and  the 
angular  constant  array.  MULT  multiplies  a  square  matrix 
times  a  column  matrix,  MULT  is  used  to  do  the  matrix 
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i'll! 'll  iplieat  ion  in  Equa  l  i  ons  (1 . 1 . G , 8 , 9 ) . 

Mil’ll  1  muI  AMPH2  are  used  to  debug  tho  routines.  AHPR1 
prints  out  ;i  column  matrix  and  ARPK2  a  square  matrix.  They 
are  no L  called  by  any  routine,  but  are  left  for  the  user  to 
use  if  he  decides  to  make  any  modifications. 

Cone,  .Fusion s 

This  chapter  lias  presented  the  main  routines  and  all  of 
its  subroutines  with  each  routine  being  explained  to  show 
how  the  theory  of  the  last  chapter  was  utilized.  It  was 
shown  how  the  routines  are  able  to  analyze  both  annular  and 
non-annular  wave-fronts.  Appendix  D  contains  detailed  flow¬ 
charts  of  all  of  the  routines,  and  Appendix  F  contains  a 
listing  of  the  entire  program.  Tho  software  has  been  writ¬ 
ten  such  that  is  does  not  rely  on  any  system  subroutines 
other  than  a  square  root  routine,  which  reduces  the  possibil 
ity  that  a  special  routine  is  not  available.  The  next  chap¬ 
ter  deals  with  the  validation  and  verification  of  the  rou¬ 
tines  presented  in  this  chapter. 


72 


_ 


V. 


Va I i da l ion 

Tlu'  previous  cha  p  t,  o  i';  ■  prosen  1  oil  (Jio  reasons  for  devel¬ 
oping;  a.  now  so  1'i.ivnre  pacliago,  the  theory  behind  the  develop 
im'ii  t  ,  and  the  actual  soriu'a  re,  This  chapter  deals  with  l  he 
validation  o  I  the  soil.  ware.  The  validation  process  is  lire!: 
en  up  into  three  sections.  The  first  section  tests  the 
orthogonal  i  ty  of  the  new  polynomials  as  the  obscuration 
ratio  changes.  The  next  section  verifies  the  software's 
ability  to  analyse  non-annular  wave-fronts.  The  final  sec¬ 
tion  deals  with  annular  wave-fronts.  Thus  all  of  the  as¬ 
pects  of  the  software  will  have  been  tested  and  verified. 

Or  thorona  1  i  t  y  Ve  r  if  ion  i  i_<  >  n 

The  major  aspect  of  this  thesis  :is  the  generation  of 
an  orthogonal  set  (>f  pol ynom ia 1  •  .  The  polynomials  must  be 
orthogonal  over  an  annular  wave-front.  This  section  veri¬ 
fies  that  such  a  set  lias  boon  genera  ted.  Since  both  the 
Zorn  ike  and  new  polynomials  use  the  same  angular  functions, 
only  t  ho  radial  portions  are  tested,  as  shown  in  Equations 
(3. 2.-1)  an  I  (3.1.23). 

The  actual  test  consisted  of  performing  the  integrals 
in  Equations  (3.2.1)  and  (3.d.22)  with  the  computer.  The 
integrals  were  approximated  by  using  Simpson's  rule.  The 
region  o  t  .integral  ion  was  divided  im.o  2,000  steps  for  the 
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initial  i  nteg  r;i  t  i  on  .  This  was  doubled  until  the  difference 
between  iterations  was  less  than  0.01.  To  show  a  compari¬ 
son  between  the  Zorn ike  and  new  radial  polynomials,  the 
first  .integration  had  limits  of  0.0  to  0.90000099,  The 
upper  limit  was  not  1.0  because  the  routine  GAMSUB  would 
try  to  compute  some  indefinite  terms. 

The  test  integrated  all  the  possible  combinations  of 
Zorn ike  radial  polynomials  and  all  the  possible  combination 
of  the  new  radial  polynomials.  After  each  integration,  the 
lower  limit  would  increase  by  0 . 05  until  the  limits  of  in¬ 
tegration  were  0.05  to  0.90099999.  Figure  19  is  a  plot  of 
the  non-zero  terms  of  the  orthogonality  condition.  All  of 
the  terms  where  the  Kronoekor  delta  is  supposed  to  be  zero 
oscillated  around  zero.  The  value  of  these  terms  were  very 
dependent  on  tin;  number  of  steps  in  the  integration  approx¬ 
imation;  therefore,  only  the  cases  where  the  Kronecker 
della  is  one  were  plotted .  Figure1  19  clearly  shows  that 
the  values  of  the  new  radial  polynomials  remained  constant 
regardless  of  the  integration  limits. 

To  verify  the  cases  where  the  Kronecker  delta  is  zero, 
the  .integrations  were  done  by  hand  using  a  variable  as  the 
lower  limit.  All  of  the  integrals  verified  the  orthogonal¬ 
ity  condition. 

This  test  verified  the  inner  product  coefficients, 
the  program  GAMBUB,  and  the  theory  presented  in  Chapter  III 


With  those  results,  t.iio  anal  ysi  s  of  wave-  fronts  can  be 
Yi.c  i  f  Led . 


Non -annul  ar  Ve i-ificat  Lon 

To  test  the  softwato's  ability  to  analyze  wave-fronts, 
two  types  of  wave- (Vents  were  used.  One  set  consisted  of 
generated  wave-fronts,  and  the  other  ones  wore  actual  wave- 
fronts  supplied  by  AFWl..  (loth  sets  were  non -annular  wave- 
fronts.  By  using  non -annular  wave-fronts,  the  program's 
ability  to  just  analyze  wave-fronts  was  tested.  With  non- 
annular  wave-fronts,  the  obscuration  ratio  is  zero  which 
results  in  the  two  sets  of  polynomials  being  identical. 

The  generated  wave-fronts  came  from  a  program  which 
had  input  parameters  consisting  of  the  number  of  coeffi¬ 
cients,  their  values,  the  center  coordinates,  and  the 
radius.  These  wave-fronts  verified  the  operation  of  the 
subroutines  FAPLR,  and  VALID.  These  subroutines  are  the 
ones  which  determined  the  type,  center,  and  radius  of  the 
wave-front.  With  the  ability  to  set  the  coefficients  to 
specific  values,  the  orthogonality  of  the  polynomials  could 
also  be  tested.  Figure  3  in  Chapter  II.  was  generated  us¬ 
ing  this  routine. 

In  every  case,  FAPEll  and  VALID  were  able  to  find  the 
center  and  radius  of  the  wave-fronts.  This  was  due  to  the 
uniformity  of  the  generated  wave-fronts.  Figure  20  shows 
the  output  when  the  second  coefficient  of  the  generated 
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Figure  20.  Output  of  Software  with  Generated  Ware-front 
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wave- front  way  set.  at  0.5,  the  radius  at  14.0,  and  the  cen¬ 


ter  at  16.5,10.5,  Figure  21  is  the  results  of  solving  for 
11  and  22  coefficients.  These  two  figures  show  the  standard 
output  of  Lho  program  written  for  this  thesis. 

As  the  number  of  coefficients  increased,  the  RMS  error 
increased  from  the  introduction  of  extraneous  coefficients. 
These  coefficients  come  about  by  the  approximation  of  contin¬ 
uous  functions  over  a  discrete  region.  As  shown  in  Figures 
20  and  21,  the  RMS  error  changes  from  0.000033  to  0.027794 
waves.  This  four  order  of  magnitude  difference  appears  to 
be  significant  error;  however,  dividing  by  the  wave  length, 
the  RMS  error  ranges  from  0.00031%  to  0.26%.  Thus  even 
though  there  is  a  large  change,  the  deviation  from  the  true 
wave-front  is  small.  Figures  22-26  show  the  results  of  the 
analysis  of  a  generated  wave-front.  The  RMS  error  for 
these  plots  was  0.026  waves.  The  difference  plots  in  Fig¬ 
ures  25  and  26  are  the  difference  between  the  actual  wave- 
front  and  the  respective  estimated  wave-front. 

After  running  over  30  tost  cases,  the  RMS  error  never 
exceeded  0.050  waves  when  preset  coefficients  were  within 
the  range  of  the  computed  coefficients.  For  example  when 
only  the  eleventh  coefficient  was  set  at  0.25,  the  RMS  er¬ 
ror  was  0.24  waves  when  the  first  six  coefficients  wore 
found.  The  RMS  error  dropped  to  0.024  waves  when  solving 
for  11  coefficients.  Thus  with  generated  non-annular 
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Igure  22,  Test  of  Non-armular  Wave-front  Analysis 
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Figure  24.  Test  of  Kon-annular  Wave-front  Analysis 


TT 


wave-fronts,  the  softWMVO  i  s  capable  of  accurately  finding 
the  coefficients  wi  th  a  minimal  armnount  of  error. 

The  second  lost  of  tin;  software  with  non-annular  wave- 
fronts  involved  data  provided  by  AFV.'L.  The  data  consists 
of  20  sequential  frames  of  wave- front.  The  first  frame  is 
circular  and  centered  in  the  colled  ion  array .  As  time 
progressed  a  knife-edge  was  inserted  in  the  wave-front 
blocking;  off  the  positive  X  data  points.  The  knife-edge 
entered  at  the  Lop>  of  quadrant  I  and  moved  in  a  -Y  direction 
until  the  wave-front  became  a  semi-circle  at  about  the 
twelfth  frame  of  data.  The  knife-edge  was  then  extracted 
in  the  opposite  direction  until  by  the  twentieth  frame  the 
wave-front,  is  nearly  circular.  The  RMS  error  from  six  co¬ 
efficients  is  plotted  in  Figure  27.  This  plot  also  has  the 
percentage  of  the  wave-front  which  was  obscured  by  the  knife- 
edge.  Tills  figure  shows  the  correspondence  between  the  RMS 
error  and  the  ammount  of  wave-front  present.  In  this  par¬ 
ticular  set  of  data,  the  wave-front  lias  a  large  ammount  of 
45°  Astigmatism  (Zorn ike  polynomial  6).  As  the  +X  data  was 
obscured ,  the  coefficient  got  larger.  It  was  this  coeffi¬ 
cient  which  caused  t  lie  RMS  error  to  .increase.  Figures  28 
to  33  show  the  frame  being  analyzed,  the  estimated  wave- 
front,  and  the  difference  between  the  two.  Figures  28  to 
30  correspond  to  frame  one  in  Figure  27  and  Figures  31  to 
33  correspond  to  frame  twelve  in  Figure  27.  Only  the 
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Figure  32.  Frame  12,  Non-annular  Wave-front  Estimate 


Zernike  arc  shown  since  the  two  polynomials  are  the  same. 
The  peaks  around  the  perimeter  of  the  difference  plots 
come  about  by  the  fitting  of  circular  regions  to  arbitrary 
regions.  The  peaks  are  not  included  in  the  RMS  error  cal¬ 
culations. 

This  section  has  dealt  with  the  analysis  of  non-annu- 
lar  wave-front.  It  has  been  shown  that  the  accuracy  of  the 
results  is  very  dependent  on  the  percentage  of  wave-front 
the  software  is  able  to  utilize. 

Annular  Verification 

This  section  deals  with  the  programs  ability  to  ana¬ 
lyze  annular  wave-fronts.  This  area  will  again  be  split 
into  two  types  of  data:  generated  and  actual.  In  the  case 
of  generated  wave-fronts,  the  same  routine  is  used  which 
generated  the  non-annular  wave-fronts.  The  only  difference 
is  the  obscuration  ratio  is  an  input  parameter;  therefore, 
any  size  of  annular  wave-front  can  be  generated.  Using 
these  generated  wave-fronts,  the  program  was  able  to  find 
and  determine  the  basic  parameters  of  the  wave-front  in  all 
cases,  even  wave-fronts  which  were  not  centered.  The  only 
time  the  software  made  a  mistake  was  when  the  wave-front 
filled  only  one  quarter  of  the  array  and  was  not  centered. 
In  this  position,  FAPER  found  the  outside  edge  first  and 
assumed  the  wave-front  was  non-annular.  Since  the  oper¬ 
ator  can  correct  problems  like  this,  the  parameters  were 


changed  and  the  program  was  able  to  find  the  coefficients 
accurately. 

The  software  was  able  to  find  the  values  of  the  co¬ 
efficients  with  the  same  degree  of  accuracy  as  that  of 
non-annular  wave-fronts.  The  closer  the  obscuration  ratio 
was  to  the  actual  value,  the  closer  the  RMS  error  of  the 
two  sets  of  polynomials.  For  example,  when  a  non-annular 
wave-front  was  analyzed  as  an  annular  one,  the  RMS  error  of 
the  new  polynomials  were  about  twice  that  of  the  Zernike 
polynomials.  As  the  obscuration  ratio  decreased,  the  RMS 
error  of  the  new  polynomials  would  approach  the  RMS  error 
of  the  Zernike  polynomials.  Thus  the  best  results  were  ob¬ 
tained  when  the  actual  and  estimated  obscuration  ratios 
matched.  The  same  results  were  seen  when  the  input  wave- 
front  was  annular.  The  RMS  error  of  the  new  polynomials; 
however,  when  the  obscuration  ratio  exceeded  about  0.6  this 
was  not  always  the  case.  In  this  regime  the  program  had 
difficulty  computing  the  coefficients  of  the  polynomials. 
Just  as  the  RMS  error  increased  in  the  non-annular  wave- 
front  when  it  was  obscured  with  a  knife-edge,  the  same 
happened  with  the  RMS  error  as  the  obscuration  ratio  in¬ 
creased  above  0.6. 

It  must  be  explained  that  the  Zernike  coefficients  are 
not  obtained  from  analyzing  the  wave-front,  but  from  con¬ 
verting  the  new  polynomial  coefficients  to  Zernike  coeffi¬ 
cients.  If  the  Zernike  coefficients  were  used,  the 


wave-front  would  always  have  to  be  assumed  to  bo  non-an- 
nular.  If  the  Zernike  polynomials  were  used  on  an  annular 
wave-front,  their  lack  of  orthogonality  would  result  in 
severe  cross-coupling.  Therefore  even  though  either  set 
may  give  a  slightly  better  result  in  the  RMS  error,  the 
main  point  is  the  wave-front  was  analyzed  with  an  orthog¬ 
onal  set  of  polynomials. 

When  actual  annular  wave-fronts  from  AFWL  were  ana- 
alyzed,  the  biggest  problem  was  entering  the  proper  par¬ 
ameters.  The  wave-front  used  was  off-center  and  part  of 
the  beam  was  obscured.  Figure  34  shows  the  annular  wave- 
front.  This  was  very  difficult  to  get  good  results  from, 
since  the  software  would  only  use  a  small  portion  of  the 
wave-front.  The  best  results  were  obtained  when  the  out¬ 
side  radius  was  increased  to  include  all  of  the  points  of 
the  array.  When  the  radius  was  increased  from  9.25  to  17.0 
and  the  obscuration  ratio  was  changed  accordingly,  the  RMS 
error  dropped  from  2.88  to  1.07  waves.  This  error  was  fur¬ 
ther  reduced  when  the  center  was  moved  to  the  approximate 
center  of  the  data,  and  the  obscuration  ratio  was  dropped 
from  0.67  to  0.15  waves.  This  change  resulted  in  the  RMS 
error  dropping  to  0.31  waves  or  2.92%. 

Thus  the  software  is  able  to  analyze  annular  wave- 
fronts  as  designed.  It  is  able  to  find  the  coefficients 
of  the  new  polynomials  and  convert  them  to  the  correspond¬ 
ing  Zernike  coefficients.  The  next  chapter  presents  the 
recommendations  and  conclusions  of  this  thesis. 


VI. 


Conclusions  and  Recommendations 


This  thesis  has  presented  a  method  of  analysing  annu¬ 
lar  wave-fronts.  Starting  with  a  brief  description  of  the 
problem  in  Chapter  I,  the  second  chapter  presented  some  in¬ 
sight  into  how  a  LWA  works.  Next  the  Patterson  and  Zernike 
polynomials  wore  compared  and  their  respective  faults  were 
presented.  Chapter  III  gave  the  reader  the  theory  behind 
the  construction  of  an  orthogonal  set  of  polynomials,  fin¬ 
ishing  with  a  means  of  converting  Zernike  polynomials  to  a 
new  set  of  orthogonal  polynomials  and  vice  versa.  Chapter 
IV  presented  the  development  of  the  software  necessary  to 
analyze  annular  wave-fronts,  and  Chapter  V  verified  the 
software's  operation. 

In  summary,  the  LWA  software  package  developed  here  is 
able  to  find  the  coefficients  of  up  to  and  including  the 
first  22  Zernike  polynomials.  The  wave-front  being  analyzed 
can  range  in  shape  from  circular  to  an  annulus  with  an  ob¬ 
scuration  ratio  of  over  0.95.  The  best  results  were  ob¬ 
tained  when  the  obscuration  ratio  did  not  exceed  0.60.  The 
software  performed  best  when  the  true  center,  obscuration 
ratio,  and  outside  radius  were  known.  With  these  parameters 
coupled  with  the  total  wave-front  area  present,  the  result¬ 
ing  RMS  error  was  consistently  loss  than  0.040  waves. 
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This  study  has  met  all  of  its  goals  which  wore: 

(1)  using  an  orthogonal  set  of  pol ynomials ,  find  the  coeffi¬ 
cients  describing  the  wave- front ,  (2)  these  coefficients 
must  not  change  as  the  number  of  coefficients  changes, 

(3)  find  the  coefficients  as  fast  as  possible.  On  a  CDC 
6600  the  software  was,  on  the  average,  able  to  find  6, 
then  11,  then  22  coefficients  of  one  frame  of  data  in  20 
seconds  of  CPU  time.  With  these  results  in  mind,  the  next 
section  considers  possible  improvements  to  the  software. 

Re  common  d  a  t i on s 

The  current  software  assumes  the  wave-front  is  either 
circular  or  annular.  One  way  to  improve  the  software  would 
be  to  have  the  software  define  the  region  and  then  generate 
a  set  of  polynomials  which  are  orthogonal  over  the  defined 
region.  This  would  allow  the  software  to  use  all  of  the 
valid  data  points.  This  could  be  further  modified  by  weight¬ 
ing  the  phase  data  values  by  the  intensity  at  each  point. 

In  so  doing,  those  phase  values  whose  corresponding  inten¬ 
sity  is  very  small  would  have  less  importance  than  those 
with  large  intensity;  thereby  giving  less  emphasis  to  phase 
data  the  system  was  just  able  to  detect. 
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V/a vo- front  Construction  from  Delta  Phase  Data 
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Figure  35.  Simple  Wave-front  Collection  Array 

The  process  of  finding  the  wave-front  from  the  delta 
phase  measurements  involves  using  the  least-squares  tech¬ 
nique.  Figure  35  shows  a  simple  collection  array  made  up 
of  nine  points.  The  a^'s  and  Ik's  are  the  X  and  Y 
delta  phase  measurements  respectively.  The  least-squares 
method  of  finding  the  wave-front  involves  minimizing  the 
square  of  the  difference  between  the  estimated  phase  value 
at  each  point  (W^)  and  the  measured  phase  difference  (a^  or 
b^.  For  the  simple  array  in  Figure  35,  the  following 
equations  result; 

o i  =  (W2-W,-a0)7  +  ( W  4 -W i -b  q ) ~ 


(A.l) 


02  -  (W2-W i -a o ) 2  +  (W3-W2-ai)2  +  ( W 5 -W 2 -b 2 ) 2 


o3  »  (Wa-Wz-aO2  +  (W6-W3-b4): 


(A. 2) 


(A. 3) 


o4  =  (w4-Vv'i-b0)2  +  (V/7-W,,-bi)2  +  (W5-W4-a2 )2 

(A. 4) 

o5  =  (We-Wa-aj)2  +  OV5-W4-a2)2  +  (W5-W2-b2)2 


+  (W0-W5-b3) 


(A. 5) 


o6  =  (W6-W3-b4)2  +  (We-Ws-aa)2  +  (Ws-We-bj)2 

(A. 6) 


o 7  =  ( W  7 -W  4 -b i ) 2  +  (W8-W7-a4 )2 


(A. 7) 


a8  =  (W8-W7-a4)  +  (W9-W7-a5)  +(W8-Ws-b3) 

(A. 8) 


a9  =  (W<,-W8-a5)  +  (V/9-W6-b5) 


(A. 9) 


Thus  to  minimize  the  sum  of  the  square  of  the  differences, 
the  derivative  of  Equations  (A. 1-9)  js  with  respect  to  the 
phase  value  at  each  point,  or  9o\,/9W^  .  Each  derivative 
is  set  equal  to  zero.  Doing  this  to  Equations  (A. 1-9)  and 
simplifying  yields: 


2W i -W2 -W4  =  -a0-bo 
3VV  2  -W  i  -W  3  -W  5  =  a0-ai-b2 


(A. 10) 
(A. 11) 


* 


2W3 -W2 -Wo 


ai-b„ 


(A. 12) 


3Wh-Wi-Wi~W7  »  b0-b,-a2  (A. 13) 

4  W  5  -Vv  ?  -W  4  -  Wt  ~\i  8  =  a2  -  a  3  +b2  —  b  3  (A. 14) 

3W  6  -W  3  -W  5  ~W  9  =  a3+b4-b5  (A. 15) 

2W  7 -W  4  -W  0  =  b3-a4  (A. 16) 

3W8-W5-W7-W9  =  a4+b3-a5  (A. 17) 

2W9-W6~We  =  a5+b5  (A. 18) 


Putting  Equations  (A. 10-18)  into  matrix  form  yields 
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This  matrix  can  be  expressed  as  AW  =  B  To  find  the 
values  of  the  phase  at  each  point,  the  matrix  A  is  in¬ 
verted  and  multiplied  times  B  .  The  matrix  A  is 
Hermitian,  irreducibly  diagonally  dominant,  and  all  of  the 
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diagonal  entries  are  positive  real  numbers  (Ref. 5)  regard¬ 
less  of  the  order  of  the  array.  Usually  one  of  the  points 
in  the'  phase  array  is  set  to  zero  as  a  starting  point. 
Following  the  procedure  outlined  above,  the  LWA  performs 
the  same  operation  on  a  32  by  32  collection  array  to  con¬ 
struct  the  wave-front  used  by  the  software  in  this  thesis. 
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Appendix  B 

Plots  of  the  First  22  Zernike  Polynomials 


The  next  22  pages  show  the  reader  what  the  first  22 
terms  of  the  Zernike  polynomials  look  like  individually. 
As  pointed  out  in  Table  II,  the  first  ten  polynomials 
are  familiar  to  those  in  the  field  of  optics.  Each  plot 
is  the  result  of  setting  the  corresponding  coefficient  to 


a  value  of  0.1  wave. 
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Appendix  C 


Zernike  Polynomial  Development 


Appendix  C 

Zorn  .ike  Polynomial  i lovolopmont 


The  basic  theory  behind  the  development,  of  the  Zorn  ike 
polynomials  comes  from  the  works  presented  by  Born  and  Wolf 
(Ref . 1 ; Appendl x  VII).  The  basic  idea  presented  by  Born  and 
Wolf  was  to  find  a  set  of  polynomials  which  are  invariant 
with  respect  to  rotation.  In  general  one  wishes  the  fol¬ 
lowing  condition  to  be  met 


//  V’a*(  X ,  ,v  )Vj}(  X ,  y  )d>:dy  -  (C.l) 

x:!+.\  ^1 

where  V  is  the  general  polynomial,  the  denotes  com¬ 

plex  conjugate,  6  is  the  Kroneeker  delta,  and  A  is  a 
normalization  constant  (Ref. 1:707).  The  rotational  invar¬ 
iance  is  expressed  as 

V ( x , y )  -  G(0)V(x',y')  (C.2) 


where 


x'  =  x  cos  0  +  y  sin  0 
y'  -  x  sin  0  +  y  sin  0 


(C.  3) 


and  G(0)  is  continuous  with  a  period  of  2tt  and  G(0)  =  1. 
G  must  also  satisfy  the  condition  of 


G(0i)C.(O2)  =  G (  0  j  +  0 2  ) 


(C.4) 


since  the  rotation  through  one  angle,  then  the  next  angle 
is  the  same  as  the  rotation  through  the  sum  of  both  angles, 
A  simple  solution  for  a  function  to  represent  0(6)  is 


G  ( 0  )  =  e 


(C.  5  ) 


where  1  is  an  integer.  Putting  Equation  (C.5)  into  Equa¬ 
tion  ( C . 2 )  and  setting  x'=r  and  y'=0  ,  and  using  Equa¬ 

tion  (C.3)  yields 


V(r  cosO,r  sinO  )  =  U(r)eJ 


(C.6) 


where  R(r)  is  strictly  a  function  of  r  .  Horn  and  Wolf 
then  show  how  if  V  is  a  polynomial  of  degree  n  then  R 

must  have  the  same  degree  n  but  with  a  power  of  r  not 
less  than  1  .  Thus,  "the  set  of  the  Zernike  clrcJ  c  poly¬ 

nomial  s  is  distinguished  from  all  other  such  sets  by  the 
property  that  it  contains  a  polynomi al  for  each  pair  of 
permisw  j.b]  c  va  1  ues  o  f  n  ( degree )  a n d  1  (angular  dependence  )  .  " 
(Ref . 1 : 708 ) (Their  emphasis)  Thus  V  can  be  expressed  as 


VX(r  cosO,r  sinO)  =  RX(r)e'XX® 
n  n 


(C.7) 


From  Equation  (C.1,7),  R  (r)  must  satisfy  the  relation 


RT(r)RX.(r)rdr  =  i\l  & 
n  '  n  n  n  n 


(C.8) 


>  -V 


whore 


1  ^ 

1  n 

a  ==  Tj- 

ii  2  k 


With  this  result,  the  basic  form  of  1  he  Zcrnil.e  polynomial 
is  expressed. 

Born  and  Wolf  then  find  an  explicit  expression  for 


Rn ( r )  as 
n 


Rra(r)  -  tm/2Q  (t) 

n  ui-m 

2 


(C.9) 


where  t=r  and  Qn_ri(t)  is  a  polynomial  in  t  of  de- 

2 

groe  Kn~m).  Since'  Q  must  still  satisfy  the  condition 
of  R,  Q  must  satisfy  Equation  (C.8)  or 


tmQk(t)Qk.(t)dt  =  a”<5kk. 


(C.10) 


where  k=£(n-m)  and  k'*=J(n^-m)  .  Born  and  Wolf  then 

equate  Q  with  the  Jacobi  polynomials.  Through  substitu¬ 
tion  and  normalization  procedures,  the  final  form  of  the 
Radial  polynomial  is 


U^r) 


£(n-m) 

E 

s=0 


(-1)* 


(n-s )  !  r 


n-2s 


s ! ( 1 ( n+m ) -s ) ! ( i ( n-m ) -s ) ! 


(C.ll) 


Thus,  using  the  explicit  expression  for  R  in  Equation  (C. 
11)  and  expanding  the  angular  function  G  into  sine  and 
cosine  terms,  the  set  of  Zernike  polynomials  can  be  found. 
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Appendix  D 


Flowed a rt s  of  Software 

Contained  in  this  appendix  are  flowcharts  of  all  of 
the  routines  .in  the  software  used  in  the  analysis  of  an¬ 
nular  wave— fronts .  The  flowcharts  appear  in  the  follow¬ 
ing  order:  Main  program,  GAMSUB,  TAPER,  VALID,  EDGE,  ZHAD, 
ZANG,  CONTUR,  RMSERR,  SERPENT,  INVERT,  MULT,  ARPR1 ,  and 
ARPR2 .  The  explanation  of  these  routines  appears  in 
Chapter  IV,  and  the  listing  of  the  program  is  in  Appendix 
E.  Figure  36  defines  all  of  the  symbols  used  in  the  flow¬ 
charts.  The  first  appearance  of  a  variable  will  be  de¬ 
fined  in  the  upper  left  hand  corner  of  the  page  it  first 


appears  on. 
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Figure  37-2.  Flowchart  of 


Firrurc  37-3.  Plowrlnrt  of  Main 


Figure  37-4.  Flowchart  of  Main 


Figure  37-8.  Flowchart  of  Main 


new 


Figure  37-10.  Flowchart  of  KnJn 


Figure  19 -Q .  Flov.'clmrt  of  FAPL'-R 


Figure  39-9.  Flowchart  of  FAPKB 


Figure  in.  Flowchart  of  FDGE 


Figure  42,  Flowchat  of  ZRAD 


Figure  43.  Flowchart  of  ZAIIG 


Figure  44.  Flowchart  of  COFTUH 
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,  ERRORS  SORT (SUU  /  NETS) 


Re  turn 


Figure  45.  Flowchart  of  RES ERR 


XC:-X  Center  (  Star 

ycj=y  »  v - T 

RADIUS :=  Outside  Radius 

BETA:-  Obscuration  Ratio _ 

INF:=  Printer  Flag  Print 


Start  )  ( XC, IC, RADIUS, B ETA, IL'F) 


Return 


Figure  46.  Flowchart  of  SEErRNT 


lCj 


Figure  1 \? •  Flowchart  of  AP.PR1 


Cstart  3  ( A ,  IS  IZE ,  NAME ) 

(  Print  ~] 

NAME 


Figure  50.  Flowchart  of  ARPR2 


Appendix  E 

Major  Arrays  General  od  by  the  Software 
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*> 


M ;gi or  Array.- ■:  Genera t ed  by  the  Fl of  twa re 

This  appendix  contains  the  arrays  generated  by  the 
software.  The  arrays  are  shown  when  the  software  is  solv¬ 
ing  for  all  22  coefficients. 
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igure  51.  Array  Created  by  ZANG 
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